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Preface 


Audience 


This manual is for compiler writers, system programmers, and application 
programmers who do not have high-level language support for Compaq 
Portable Mathematics Library (CPML)! routines in their language of choice. 
This audience needs to access CPML routines directly from their application 
programs. 


Document Organization 


This manual consists of the following: 


Chapter 1 gives a general overview of the mathematics library and discusses 
supported data types, exception behavior, and IEEE considerations. 


Chapter 2 explains the presentation format of a CPML routine and how to 
interpret a routine’s interface, and alphabetically lists the routines. 


Appendix A lists the floating-point boundary values used by the CPML 
routines. 


Appendix B contains the complete list of entry-point names. 
The Glossary lists mathematical terms and symbolic names used in this 
manual, and provides a brief definition. 

Related Documents 


Some books in Compaq’s documentation sets help meet the needs of several 
audiences. For example, the information in some system books is also used by 
programmers. Keep this in mind when searching for information on specific 
topics. 


1 Formerly DIGITAL Portable Mathematics Library (DPML). 
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Use the documentation overview and the master index information for 
your operating system when searching for hardcopy information on a topic. 
They provide information on all of the books in your operating system's 
documentation set. 


OpenVMS Documentation 
For additional information on the Open Systems Software Group (OSSG) 
products and services, access the following OpenVMS World Wide Web address: 


http: //www.openvms.digital.com 


Compaq Tru64 UNIX Documentation 

The printed version of the Compaq Tru64 UNIX documentation set uses 
letter icons on the spines of books to help specific audiences quickly find the 
books that meet their needs. (You can order the printed documentation from 
Compaq.) The following list describes this convention: 


Audience Icon 
General users G 
System and network administrators S 
Programmers P 
Device driver writers D 
Reference page users R 


For additional information about Compag Tru64 UNIX products and services, 
access the following World Wide Web address: 


http: //www.unix.digital.com 


Reader’s Comments 


viii 


Compaq welcomes your comments on this manual. A Reader’s Comments form 
is located online in the following location: 


UNIX /usr/doc/readers_comment.txt 
OpenVMS SYS$HELP:OPENVMSDOC_COMMENTS.TXT 


Print or edit the online form and send us your comments by one of the 
following methods: 


Internet dpml@digital.com 


Fax (603) 884-0120, Attention: Math Run Time Library Documentation, 
ZK O2-3/Q18 


Mail Compaq Computer Corporation 
Math Run Time Library Documentation, ZK O2-3/Q18 
110 Spit Brook Rd. 
Nashua, NH 03062-2698 
USA 


Please include the following information along with your comments: 


¢ The full title of the book and the order number. (The order number is 
printed on the title page of this book and on its back cover.) 


e The section numbers and page numbers of the information on which you 
are commenting. 


¢« The operating system and version of Compaq software that you are using. 
For example, Tru64 UNIX Version 5.0, or OpenVMS Version 7.1. 


e If known, the type of processor that is running the Compaq software. For 
example, AlphaServer 2000. 


Compaq publications groups cannot respond to system problems or technical 
support inquiries. Please address technical questions to your local system 
vendor or to the appropriate Compag technical support office. Information 
provided with the software media explains how to send problem reports to 
Compaq. 


How To Order Additional Documentation 


Additional documentation for your operating system is available by contacting 
Compaq. For help in ordering or in deciding which documentation best meets 
your needs, call 800-344-4825. 


For OpenVMS, you can also use the following World Wide Web address to order 
documentation: 


http: //www.openvms.digital.com:81/ 


Conventions 


In this book, every use of OpenVMS means the Compag OpenVMS operating 
system, and every use of UNIX or Tru64 UNIX means the Compaq Tru64 
UNIX operating system. 


The following conventions are used in this manual: 


Ctri/x 


PF1 x 


Return 


{} 
bold text 


italic text 


A sequence such as Ctrl/x indicates that you must hold 
down the key labeled Ctrl while you press another key or a 
pointing device button. 


A sequence such as PF 1 x indicates that you must first press 
and release the key labeled PF 1 and then press and release 
another key or a pointing device button. 


In examples, a key name enclosed in a box indicates that 
you press a key on the keyboard. (In text, a key name is not 
endiosed in a box.) 


A horizontal ellipsis in examples indicates one of the 
following possibilities: 


¢ Additional optional arguments in a statement have 
been omitted. 


¢ The preceding item or items can be repeated one or 
more times. 


e Additional parameters, values, or other information can 
be entered. 


A vertical ellipsis indicates the omission of items from a 
code example or command format; the items are omitted 
because they are not important to the topic being discussed. 


In command format descriptions, parentheses indicate that, 
if you choose more than one option, you must enclose the 
choices in parentheses. 


In command format descriptions, brackets indicate optional 
elements. You can choose one, none, or all of the options. 
(Brackets are not optional, however, in the syntax of a 
directory name in an OpenVMS file specification or in 

the syntax of a substring specification in an assignment 
statement.) 


In command format descriptions, braces indicate a required 
choice of options; you must choose one of the options listed. 


This text style represents the introduction of a new term or 
the name of an argument, an attribute, or a reason. 


Italic text indicates important information, complete titles 
of manuals, or variables. Variables include information 
that varies in system output (Internal error number), in 
command lines (/PRODUCER=name), and in command 
parameters in text (where devicename contains up to five 
alphanumeric characters). 


UPPERCASE TEXT 


Monospace text 


numbers 


Uppercase text indicates a command, the name of a routine, 
the name of a file, or the abbreviation for a system privilege. 


Monospace type indicates code examples and interactive 
screen displays. In the C programming language, 
monospace type in text identifies the following elements: 
keywords, the names of independently compiled external 
functions and files, syntax summaries, and references to 
variables or identifiers introduced in an example. 


A hyphen at the end of a command format description, 
command line, or code line indicates that the command or 
statement continues on the following line 


All numbers in text are assumed to be decimal unless 
otherwise noted. Nondecimal radixes—binary, octal, or 
hexadecimal—are explicitly indicated. 


xi 
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Introduction to CPML 


The Compaq Portable Mathematics Library (referred to as CPML) includes a 
wide variety of mathematical routines that cover the following areas: 


e Floating-point trigonometric function evaluation 

e Exponentiation, logarithmic, power function evaluation 
¢ Hyperbolic function evaluation 

e Algebraic function evaluation 

* Complex function evaluation 

e« Complex exponentiation 

e Miscellaneous function evaluation 


This manual documents the CPML routines and, in particular, how they 
behave when given an exceptional input argument. It also documents 
operating system entry points and supported floating-point data types. 


1.1 Overview 


Developing software within the confines of high-level languages like Fortran 
and C greatly increases the portability and maintainability of your source code. 
Many high-level languages support mathematical function evaluation. CPML 
was developed to provide a common set of routines that support many of the 
common mathematical functions across a wide variety of operating systems, 
hardware architectures, and languages. 


In most cases, the common mathematical functions behave the same way 

for all languages and platforms. Occasionally, however, high-level language 
definitions of the same mathematical function will differ for specific input 
values. For example, in Fortran, log(-1.0) causes a program abort, while in C, 
log(-1.0) quietly returns a system-defined value. 
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This document uses the term exceptional arguments to refer to values in the 
following situations: 


e Values for which high-level languages disagree on the function behavior. 
¢« Values that are mathematically undefined or out of range. 

e Values for which the function would overflow or underflow. 

See Section 1.3 for more detail on exceptional arguments. 


To provide uniform quality of mathematical functions for all languages on 
your system, CPML traps exceptional arguments and invokes a system- 
specific routine called the CPML exception handler. The exception handler 

is designed to work in conjunction with high-level language compilers and 
run-time libraries (RTLs) to provide specific language semantics for exceptional 
arguments. This means that the user-visible behavior of a given function 
called from a given language is not necessarily determined by the routines 

in the CPML library but rather by a combination of several entities acting in 
concert. 


Note 


Compag strongly recommends that you limit your access to the CPML 
routines documented in this manual to the high-level language syntax 
of your choice, thereby guaranteeing the behavior of the routines across 
platforms. Because of the complex relationship between high-level 
languages and CPML routines, the behavior of direct calls to CPML 
routines may change from release to release. 


1.2 Data Types 


1-2 


CPML is designed to support mathematics function evaluation for multiple 
data types. These data types include integer, floating-point, and complex 
floating-point. 


The integer data type, identified as int throughout this manual, is the natural 
size signed integer for a particular platform. On a 32-bit system, int is a 32-bit 
signed integer, and on a 64-bit system, int is a 64-bit signed integer. 


The floating-point types referred to in this document are F_FLOAT, G FLOAT, 
X_FLOAT, S_FLOAT, and T_FLOAT respectively. When it is not necessary to 
distinguish between the different floating types, they are referred to collectively 
as F_TYPE. Your platform may support all or a subset of these floating-point 
data types. For example, CPML on OpenVMS Alpha systems supports the 
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following floating-point data types: VAX single- and double-precision, |EEE 
single: and double-precision, and |EEE extended-precision. CPML on Tru64 
UNIX Alpha systems supports only IEEE single and doubleprecision data 
types. Table 1-1 describes the floating-point data types. 


Table 1-1 Floating-Point Data Types 


F_TYPE Description 

S_FLOAT 32-bit IEEE single-precision number 
T_FLOAT 64-bit IEEE double-precision number 
X_FLOAT 128-bit |EEE extended-precision number 
F_FLOAT 32-bit VAX single-precision number 

G_ FLOAT 64-bit VAX double-precision number 


In addition to the data types mentioned in Table 1-1, CPML also provides 
routines that return two values of the same floating-point type; for example, 
two S_TYPE values or two G_TYPE values. In the discussion that follows, 
these pairs of floating-point data type values are referred to as F COMPLEX. 
Refer to Table 1-2. This document uses F_COMPLEX to indicate that a given 
routine returns two different values of the same floating-point data type. 


The mechanism for returning two floating-point values from CPML routines 
varies from platform to platform. However, on OpenVMS Alpha systems, 
F_ COMPLEX data is returned in consecutive floating-point registers and is 
accessible only through a high-level language, like Fortran, that specifically 
allows access to them. 


A complex number, z, is defined as an ordered pair of real numbers. The 
convention used in this manual to define an ordered pair of real numbers as 
complex is as follows: 


e The first number is the real part of the complex number. 


« The second number is preceded by i and is the imaginary part of the 
complex number. 


¢ A separator character (plus sign) is used to associate and separate the real 
and the imaginary number. 


For example: 


Z =x +iy 
z =sin x +icos y 
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CPML includes complex functions, for example, the complex sine, csin(x,y), 
defined to be sin(x +iy). Complex function routines like csin(), which have 
complex input, accept floating-point numbers in pairs and treat them as if they 
are real and imaginary parts of a complex number. 


In the previous two examples, the first floating-point values are defined by 

xX and sin x, respectively, and are the real part of the complex number. The 
second floating-point values used in the examples are defined by iy and icos 
y, respectively, and are the imaginary part of the complex number. Similarly, 
CPML routines that return complex function values return two floating-point 
values. Taken together, these two floating-point values represent a complex 
number. 


CPML supports the floating-point complex types described in Table 1-2. CPML 
complex functions can be accessed only through high-level languages that 
support the complex data type. Use only the data types supported by your 
system. 


Table 1-2 Floating-Point Complex Data Types 
F_COMPLEX Description’ 


S_FLOAT_COMPLEX An ordered pair of S FLOAT quantities, representing a 
single-precision complex number 


T_FLOAT_COMPLEX An ordered pair of T FLOAT quantities, representing a 
double-precision complex number 


X_FLOAT_COMPLEX An ordered pair of X_FLOAT quantities, representing an 
extended-precision complex number 


F_FLOAT_COMPLEX An ordered pair of F_FLOAT quantities, representing a 
single-precision complex number 


G_ FLOAT_COMPLEX An ordered pair of G FLOAT quantities, representing a 
double-precision complex number 


1The lower addressed quantity is the real part; the higher addressed quantity is the imaginary 
part. 
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Not all mathematical functions are capable of returning a meaningful result 
for all input argument values. Any argument value passed to a CPML 
routine that does not return a meaningful result, or is defined differently for 
different environments, is referred to as an exceptional argument. Exceptional 
arguments that result in an exception behavior are documented in the 
Exceptions section of each CPML routine in Chapter 2. 
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Exceptional arguments typically fall into one of two categories: 


e« Domain errors or invalid arguments. These are arguments for which a 
function is not defined. For example, the inverse sine function, asin, is 
defined only for arguments between -1 and +1 inclusive. Attempting to 
evaluate acos(-2) or acos(3) results in a domain error or invalid argument 
error. 


e Range errors. These errors occur when a mathematically valid argument 
results in a function value that exceeds the range of representable values 
for the floating-point data type. Appendix A gives the approximate 
minimum and maximum values representable for each floating-point data 
type. 


1.4 Exception Conditions and Exception Behavior 


CPML routines are designed to provide predictable and platform-consistent 
exception conditions and behavior. When an exception is triggered in a CPML 
routine, two pieces of information can be generated and made available to the 
calling program for exception handling: 


e A notification that an exception has occurred. The mechanics of exception 
notification vary from platform to platform (for example, signaling, 
trapping, set errno). 


e A return value. If your environment allows your routine to continue after 
raising an exception condition (with an exception handler for example), 
then a return value is made available upon completion of the routine. 


The exception condition-handling mechanisms on your platform dictate how 
you can recover from an exception condition, and whether you can expect to 
receive an exception notification, a return value, or both, from a CPML routine. 


The Exceptions section of each CPML routine documents each exceptional 
argument that results in an exception behavior. In addition to the exceptional 
arguments, an indication of how the CPML routines treat each argument is 
given. Exceptional arguments are sometimes presented in terms of symbolic 
constants. 


For example, the following table lists the exceptional arguments of the 
exponential routine, exp(x): 


Exceptional Argument Exception Condition/Routine Behavior 


x >In(max_float) Overflow 
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Exceptional Argument Exception Condition/Routine Behavior 


x <In(min_float) Underflow 


The exceptional arguments indicate that whenever x >In(max_float) 
or x <In(min_float), CPML recognizes an overflow or underflow condition, 
respectively. 


The symbolic constants In(max_float) and In(min_float) represent the natural 
log of the maximum and minimum representable values of the floating-point 
data type in question. The actual values of In(max_float) and In(min_float) are 
described in Appendix A. 


CPML recognizes three predefined conditions: overflow, underflow, and invalid 
argument. Table 1-3 describes the default action and return value of each 
condition. 


Table 1-3 Default Action and Return Values for Exception Conditions 


Exception Condition Default Action Return Value 
Overflow Trap HUGE_RESULT 
U nderflow Continue Quietly 0) 

Invalid argument Trap INV_RESULT 


The values HUGE RESULT and INV_RESULT are data-type dependent. 


For IEEE data types, HUGE_RESULT and |INV_RESULT are the floating-point 
encodings for Infinity and NaN respectively. 


For VAX data types, HUGE_RESULT and INV_RESULT are max_float and 0 
respectively. 


1.5 IEEE Std 754 Considerations 
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The Institute of Electrical and Electronics Engineers (IEEE) ANSI/IEEE Std 
754-1985, IEEE Standard for Binary Floating-Point Arithmetic data types 
include denormalized numbers (very close to zero). The standard supports the 
concept of “Not-a-Number” or NaN to represent indeterminate quantities, and 
uses plus infinity or minus infinity (so that they behave in arithmetic) like the 
mathematical infinities. Whenever a CPML routine produces an overflow or 
indeterminate condition, it generates an infinity or NaN value. 


All CPML routines, except one, return a NaN result when presented with a 
NaN input. The only exception is pow(NaN,0) =1 in ANSI C. 
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1.6 X/Open Portability Guide Considerations 


Table 1-4 lists the routines described in this manual that conform to the 
requirements of the X/ Open Portability Guide Version 4 (XPGA4), or are 
implemented as UNIX extensions to the XPG4 standard (XPG4-UNIX). 
Descriptions of these routines appear in Chapter 2 under the generic 
function name listed in Table 1-4. Platform-specific entry-points are listed in 
Appendix B. 


Table 1-4 XPG4 Conformant Routines 


Routine Conforms to Standard Generic Function Name 
acos XPG4 acos 
acosh XPG4-UNIX acosh 
asin XPG4 asin 
asinh XPG4-UNIX asinh 
atan XPG4 atan 
atan2 XPG4 atan 
atanh XPG4-UNIX atanh 
ceil XPG4 ceil 
cos XPG4 cos 
cosh XPG4 cosh 
cot XPG4 cot 
erf XPG4 erf 
erfc XPG4 erf 
exp XPG4 exp 
expm1 XPG4-UNIX exp 
fabs XPG4 fabs 
floor XPG4 floor 
fmod XPG4 fmod 
frexp XPG4 frexp 
gamma XPG4 Igamma 
hypot XPG4 hypot 
ilogb XPG4-UNIX ilogb 


(continued on next page) 
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Table 1-4 (Cont.) XPG4 Conformant Routines 


Routine Conforms to Standard Generic Function Name 
isnan XPG4 isnan 

jo XPG4 bessel 

jl XPG4 bessel 

jn XPG4 bessel 
Idexp XPG4 Idexp 
Igamma XPG4 Igamma 
log XPG4 log 

log10 XPG4 log 

loglp XPG4-UNIX log 

logb XPG4-UNIX logb 
modf XPG4 modf 
nextafter XPG4-UNIX nextafter 
pow XPG4 pow 
remainder XPG4-UNIX remainder 
rint XPG4-UNIX rint 
scalb XPG4-UNIX scalb 

sin XPG4 sin 

sinh XPG4 sinh 

tan XPG4 tan 

tanh XPG4 tanh 

yO XPG4 bessel 

yl XPG4 bessel 

yn XPG4 bessel 
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CPML Routines 


Each CPML routine documented in this chapter is presented in the following 
format: 


e Routine name—A brief name to identify the function of the routine A 
routine may contain more than one function. 


e |nterface—What the routine expects to receive and what it returns. See 
Section 2.1 for more information. 


¢ Description—Additional information, including the permitted range of 
input values and generic calculations used to compute the results. 


e Exceptions—A description of how the routine behaves when given a specific 
exceptional input argument. 


2.1 CPML Routine Interface 


The interface to each function is: 
RETURN_TYPE generic_interface_name (INPUT_ARG_TYPE...) 


Each of these is described below. 


RETURN_TYPE 

The data type of the value returned by the routine to your application program. 
Each routine returns a specific class of data type. For example, either F_TYPE 
or F_ COMPLEX can appear in a CPML interface described in Chapter 2. The 
supported data types are described in Section 1.2. 


generic_interface_name 

The generic name. CPML routines in this chapter are listed in alphabetic 
order by their interface names. Some CPML routines may be available in the 
syntax of your high-level language. Fortran and C are examples. To maximize 
the portability of your application, use the corresponding mathematical routine 
described in your high-level language, and directly call only the routines 
documented in this manual that are not supported by your language. Refer to 
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Appendix B for the specific entry-point names needed to directly call a CPML 
routine from your platform. 


INPUT_ARG_TYPE... 

The number and type of input arguments provided by your application. Some 
routines require more than one argument. Arguments must be coded in the 
order shown in the interface section of each routine described in this chapter. 
The supported data types for arguments are described in Section 1.2. 


Note 


Unless otherwise noted, arguments are read-only and passed by value. 
Arguments passed by another mechanism are prefaced by an asterisk 
(*); for example, *n in the frexp() routine. 


2.2 Specific Entry-Point Names 


Each generic interface name documented in the interface section of a routine 
description corresponds to one or more specific entry-point names described 
in Appendix B. For example, on OpenVMS Alpha systems, the acosd function 
has five entry-point names; one for each available floating-point data type. 
The acosd entry-point names are math$acosd_f, math$acosd_s, math$acosd_x, 
math$acosd_g, and math$acosd_t. On Tru64 UNIX Alpha systems, the acosd 
function has two entry-point names corresponding to their supported data 
types: S FLOAT and T_FLOAT. The two entry-point names are acosdf for 

S FLOAT input arguments and acosd for T FLOAT arguments. Use the 
specific entry-point name that corresponds to the input argument data type. 


2.3 Working with Exception Conditions 


Each CPML routine description contains a table of exceptions. Each exception 
listed in the table represents an exceptional case that is handled in a platform- 
specific manner. For example, the atan2() exception table contains the 
following two entries: 


Exceptional Argument Routine Behavior 
y=x=0 Invalid argument 
| y| =| x| =infinity Invalid argument 


The first entry describes an exception condition containing two input 
arguments with zero values. Upon detecting this error, the routine behavior 
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signals the “invalid argument” condition. The second entry is applicable only to 
platforms supporting signed or unsigned infinity values. Here, if the absolute 
value of both input arguments is equal to infinity, an “invalid argument” 
condition is signaled. 


The exact behavior of a routine that detects an exceptional argument varies 
from platform to platform and is sometimes dependent on the environment in 
which it is called. The behavior you see depends on the platform and language 
used. It also depends on how the routine was called and the interaction of the 
various layers of software through which the call to the routine was made. 
Remember, access to a CPML routine can be made through direct access (a 
CALL statement written by a programmer in a source code statement) or 
through indirect access (from compiler-implemented mathematical syntax). 


The default behavior for detecting the x=y=0 arguments is to generate an 
exception trap when accessing atan2() indirectly through Fortran compiler 
syntax. C compiler syntax for the atan2() routine sets errno and returns a 
NaN when give the same input. In these cases, your compiler documentation 
provides you with information on how to work with exception conditions. 


2.4 CPML Routine Interface Examples 


This section discusses the atan2() and cdiv() interfaces and explains how to 
interpret them. The explanations given in this section apply to all CPML 
routines. 


2.4.1 atan2() Interface 
The interface to the atan2() routine is: 


F_TYPE atan2 (F_TYPE y, F_TYPE x) 


The routine name atan2() is the high-level language source-level name that 
gets mapped to a specific entry-point name documented in Appendix B. This 
is the name that appears in compiler documentation for this mathematical 
routine. The appropriate entry-point name is automatically selected when 
atan2() is called from high-level language syntax. This selection depends 
upon the data type of the input arguments. If you make direct calls to this 
routine, you must manually select the proper entry-point name documented in 
Appendix B for the data type of your input arguments. 


The format of the atan2() routine shows that it expects to receive two input 
arguments by value. Both arguments must be the same F_TYPE. The returned 
value will also be the same F_TYPE as the input arguments. 
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For example, on OpenVMS Alpha systems, the G_FLOAT entry-point name is 
math$atan2_g(). It takes two G_FLOAT arguments by value and returns a 
G_ FLOAT result. 


For Tru64 UNIX Alpha systems, the S FLOAT entry-point name is atan2f(). 
The routine takes two S FLOAT input arguments by value and returns an 
S FLOAT result. 


2.4.2 cdiv() Interface 
The interface to the cdiv() routine is: 


F_COMPLEX cdiv (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE qd) 


The routine name cdiv() is the generic name that gets mapped to a specific 
entry-point name documented in Appendix B. Selection of the appropriate 
entry-point name is done automatically when cdiv() is called from high-level 
language syntax. This selection depends upon the data type of the input 
arguments. Again, if you make direct calls to this routine, you must manually 
select the proper entry-point name documented in Appendix B for the data 
type of your input arguments. 


The format of the cdiv() routine shows that it expects to receive four input 
arguments by value. All arguments must be the same F_TYPE. The returned 
value will be an F_COMPLEX data type and will be the same base data type 
as the input arguments. 


For example, on OpenVMS Alpha systems, the F_FLOAT entry-point name is 
math$cdiv_f(). This routine takes four F_FLOAT input arguments by value 
and returns an F_FLOAT COMPLEX result in an ordered pair of F_ FLOAT 
quantities. 


For Tru64 UNIX Alpha systems, the S FLOAT entry-point name is cdivf(). 
This routine takes four S FLOAT input arguments by value and returns an 
S FLOAT_COMPLEX result. 
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acos - Arc Cosine of Angle 


Interface 
F_TYPE acos (F_TYPE x) 
F_TYPE acosd (F_TYPE x) 


Description 


acos( ) computes the principal value of the arc cosine of x in the interval [0,pi] 
radians for x in the interval [-1,1]. 


acosd() computes the principal value of the arc cosine of x in the interval 
[0,180] degrees for x in the interval [-1,1]. 


Exceptions 
Exceptional Argument Routine Behavior 
| x| >1 Invalid argument 
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acosh - Hyperbolic Arc Cosine of Angle 


Interface 
F_TYPE acosh (F_TYPE x) 


Description 


acosh() returns the hyperbolic arc cosine of x for x in the interval [1,+infinity], 
where acosh(x) =In(x +sqrt(x**2 - 1)). acosh() is the inverse function of cosh(), 
where acosh(cosh (x)) =x. 


Exceptions 
Exceptional Argument Routine Behavior 
x<1 Invalid argument 
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asin - Arc Sine of Angle 


Interface 
F_TYPE asin (F_TYPE x) 
F_TYPE asind (F_TYPE x) 


Description 


asin() computes the principal value of the arc sine of x in the interval 
[-pi/2,pi/2] radians for x in the interval [-1,1]. 


asind() computes the principal value of the arc sine of x in the interval [-90,90] 
degrees for x in the interval [-1,1]. 


Exceptions 
Exceptional Argument Routine Behavior 
| x| >1 Invalid argument 
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asinh - Hyperbolic Arc Sine of Angle 


Interface 
F_TYPE asinh (F_TYPE x) 


Description 


asinh() returns the hyperbolic arc sine of x for x in the interval 
[-infinity, +infinity], where asinh(x) =In(x + sqrt(x**2 +1)). asinh() is the 
inverse function of sinh(), where asinh(sinh (x)) =x. 


Exceptions 


None. 
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atan - Arc Tangent of Angle with One Argument 


Interface 
F_TYPE atan (F_TYPE x) 
F_TYPE atand (F_TYPE x) 


Description 


atan() computes the principal value of the arc tangent of x in the interval 
[-pi/2,pi/2] radians for x in the interval [-infinity, +infinity]. 


atand() computes the principal value of the arc tangent of x in the interval 
[-90,90] degrees for x in the interval [-infinity, +infinity]. 


Exceptions 


None. 
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atan2 - Arc Tangent of Angle with Two Arguments 


Interface 
F_TYPE atan2 (F_TYPE y, F_TYPE x) 
F_TYPE atand2 (F_TYPE y, F_TYPE x) 


Description 


atan2() computes the angle in the interval [-pi,pi] whose arc tangent is y/x 
radians for x and y in the interval [-infinity, tinfinity]. The sign of atan2() 
is the same as the sign of y. The atan2(y, x) function is computed as follows 
where f is the number of fraction bits associated with the data type: 


Value of Input Arguments Angle Returned 
x=0or gee2t n/2% (St gy) 
z>0Oand y/x < gftl atan. (y/=) 

a<0Oand y/x < gftl T * (St gry) + atan (y/=) 


atand2() computes the angle in the interval [-180,180] whose arc tangent is y/x 
degrees for x and y in the interval [-infinity, tinfinity]. The sign of atand2() is 
the same as the sign of y. 


Exceptions 
Exceptional Argument Routine Behavior 
y =x =0 Invalid argument 
| y| =infinity and | x| =infinity Invalid argument 
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atanh - Hyperbolic Arc Tangent of Angle 


Interface 
F_TYPE atanh (F_TYPE x) 


Description 


atanh() returns the hyperbolic arc tangent of x for x in the interval (-1,1). 
atanh() is the inverse function of tanh(), where atanh(tanh (x)) =x. 


atanh(x) is computed as 1/2 In((1+x)/(1-x)). 


Exceptions 
Exceptional Argument Routine Behavior 
|x| >or=1 Invalid argument 
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bessel - Bessel Functions 


Interface 
F_TYPE jO (F_TYPE x) 
F_TYPE ji (F_TYPE x) 
F_TYPE jn (int n, F_TYPE x) 
F_TYPE yO (F_TYPE x) 
F_TYPE y1 (F_TYPE x) 
F_TYPE yn (int n, F_TYPE x) 


Description 


jO() and j1() return the value of the Bessel function of the first kind of orders 0 
and 1 respectively. 


jn() returns the value of the Bessel function of the first kind of order n. 


yO() and y1() return the value Bessel function of the second kind of orders 0 
and 1 respectively. 


yn() returns the value of the Bessel function of the second kind of order n. 


The value of x must be positive for the y family of Bessel functions. The value 
of n specifies some integer value. 


Exceptions 
Exceptional Argument Routine Behavior 
(yO(), y1(), yn()) x <0 Invalid argument 
(yO(), y1(), yn()) x =0 Overflow 


The j1() and jn() functions can result in an underflow as x gets small. The 
largest value of x for which this occurs is a function of n. 


The y1() and yn() functions can result in an overflow as x gets small. The 
largest value of x for which this occurs is a function of n. 
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cabs - Complex Absolute Value 


Interface 
F_TYPE cabs (F_TYPE x, F_TYPE y) 


Description 


cabs(x,y) is defined as the square root of (x**2 + y**2) and returns the same 
value as hypot(x,y). 


Exceptions 
Exceptional Argument Routine Behavior 
sqrt(x**2 + y**2) >max_float Overflow 

See Also 


Appendix A, Critical Floating-Point Values 
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cbrt - Cube Root 


Interface 
F_TYPE cbrt (F_TYPE x) 


Description 


cort() returns the rounded cube root of x. 


Exceptions 


None. 
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ccos - Cosine of Angle of a Complex Number 


Interface 
F_COMPLEX ccos (F_TYPE x, F_TYPE y) 


Description 
ccos() returns the cosine of a complex number, x +iy. 


ccos(x,y) is defined as cos (x +iy) =(cos x cosh y-isin x sinh y). 


Exceptions 
Exceptional Argument Routine Behavior 
| y| =infinity Invalid argument 
(sin x sinh y) > max_float Overflow 
(cos x cosh y) >max_float Overflow 
See Also 


Appendix A, Critical Floating-Point Values 
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cdiv - Complex Division 


Interface 
F_COMPLEX cdiv (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 


Description 


cdiv() returns the quotient of two complex numbers: (a +ib)/(c +id). 


Exceptions 
Exceptional Argument Routine Behavior 
c=d=0 Invalid argument 
a=b=c=d=0 Invalid argument 
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ceil - Ceiling 


Interface 
F_TYPE ceil (F_TYPE x) 


Description 


ceil() returns the smallest floating-point integer value greater than or equal to 
x. 


Exceptions 


None. 


CPML Routines CPML-17 


CPML Routines 
cexp - Complex Exponential 


cexp - Complex Exponential 


Interface 
F_COMPLEX cexp (F_TYPE x, F_TYPE y) 


Description 
cexp() returns the exponential of a complex number. 
cexp(x,y) is defined as e**(x +iy) =e**x cos y +ie**x sin y. 


Exceptions 
Exceptional Argument Routine Behavior 
| y| =infinity Invalid argument 
| e**x cos y|_ >max_float Overflow 
| ex sin y| >max_float Overflow 
| e**x cos y| <min_float Underflow 
| ex sin y| <min_float Underflow 
See Also 


Appendix A, Critical Floating-Point Values 
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clog - Complex Natural Logarithm 


Interface 
F_COMPLEX clog (F_TYPE x, F_TYPE y) 


Description 
clog() returns the natural logarithm of a complex number. 
clog(x,y) is defined as In(x +iy) =1/2 In(x**2 + y**2) +iatan2 (y,x). 


Exceptions 
Exceptional Argument Routine Behavior 
y=x=0 Invalid argument 
| y| =| x| Snfinity Invalid argument 
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cmul - Complex Multiplication 


Interface 
F_COMPLEX cmul (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 


Description 
cmul() returns the product of two complex numbers. 
cmul(a,b,c,d) is defined as (a +ib) * (c +id). 


Exceptions 


None. 
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copysign - Copy Sign 


Interface 
F_TYPE copysign (F_TYPE x, F_TYPE y) 


Description 


copysign() returns x with the same sign as y. IEEE Std 754 requires 
copysign(x,NaN) = +x or -x. 


Exceptions 


None. 


CPML Routines CPML-21 


CPML Routines 
cos - Cosine of Angle 


cos - Cosine of Angle 


Interface 
F_TYPE cos (F_TYPE x) 
F_TYPE cosd (F_TYPE x) 


Description 
cos() computes the cosine of x, measured in radians. 


cosd() computes the cosine of x, measured in degrees. 


Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
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cosh - Hyperbolic Cosine of Angle 


Interface 
F_TYPE cosh (F_TYPE x) 


Description 
cosh() computes the hyperbolic cosine of x. 
cosh(x) is defined as (e**x + e**(-x))/2. 


Exceptions 
Exceptional Argument Routine Behavior 
| x] >In(2 * max_float) Overflow 

See Also 
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cot - Cotangent of Angle 


Interface 
F_TYPE cot (F_TYPE x) 
F_TYPE cotd (F_TYPE x) 


Description 
cot() computes the cotangent of x, measured in radians. 


cotd() computes the cotangent of x, measured in degrees. 


Exceptions 
Exceptional Argument Routine Behavior 
(cot) x=0 Overflow 
(cotd) | x] =multiples of 180 degrees Overflow 
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cpow - Complex Power 


Interface 
F_COMPLEX cpow (F_TYPE a, F_TYPE b, F_TYPE c, F_TYPE d) 


Description 
cpow() raises a complex base (a +ib) to a complex exponent (c +id). 
cpow(a,b,c,d) is defined as e**((c +id) In(a +ib)). 


Exceptions 
Exceptional Argument Routine Behavior 
sqrt (a**2 + b**2) >max_float Overflow 
c/2 * In(a**2 + b**2) >max_float Overflow 
c/2 * In(a**2 + b**2) - (d * atan2(b,c)) >max_float Overflow 
See Also 
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csin - Sine of Angle of a Complex Number 


Interface 
F_COMPLEX csin (F_TYPE x, F_TYPE y) 


Description 
csin() computes the sine of a complex number, x +/iy. 


csin(x,y) is defined as sin (x +iy) =(sin x cosh iy +icos x sinh iy). 


Exceptions 
Exceptional Argument Routine Behavior 
| y| =infinity Invalid argument 
(sinh x sin y) >max_float Overflow 
(cosh x cos y) >max_float Overflow 
See Also 
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csqrt - Complex Square Root 


Interface 
F_COMPLEX csqrt (F_TYPE x, F_TYPE y) 


Description 


csqrt() computes the square root of a complex number, x +iy. The real part of 
csqrt is greater than or equal to zero. 


csqrt(x,y) is defined as the square root of (x +iy). 


Exceptions 
None. 
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cvt_ftof - Convert Between Supported Floating-Point Data 


Types 


Interface 


int cvt_ftof void *x, int x_type, void “y, int y_type, options 


Description 


Note 


This routine does not apply to OpenVMS Alpha. OpenVMS Alpha users 
should use the CVT$FTOF routine documented in the OpenVMS RTL 


Library (LIB$) Manual. 


cvt_ftof() converts a floating-point value from one data type to another. x 
points to the input value to be converted, and y points to the converted result. 
The conversion is subject to the options specified in the options (bit field) 


argument. 


x_type and y_type identify the data type of x and y as follows: 


Values for x_type and y_type 


Floating-Point Data Type 


CVT_VAX_F 

CVT_VAX_D 

CVT_VAX_G 

CVT_VAX_H 

CVT_IEEE_S 

CVT_IEEE_T 

CVT_IEEE_X 
CVT_BIG_ENDIAN_IEEE_S 
CVT_BIG_ENDIAN_IEEE_T 
CVT_BIG_ENDIAN_IEEE_X 
CVT_IBM_SHORT 
CVT_IBM_LONG 
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VAX F Floating ( 4 bytes) 

VAX D Floating ( 8 bytes) 

VAX G Floating ( 8 bytes) 

VAX H Floating (16 bytes) 

IEEE Little Endian S Floating ( 4 bytes) 
IEEE Little Endian T Floating ( 8 bytes) 
IEEE Little Endian X Floating (16 bytes) 
IEEE Big Endian S Floating ( 4 bytes) 
IEEE Big Endian T Floating ( 8 bytes) 
IEEE Big Endian X Floating (16 bytes) 
IBM_Short_Floating ( 4 bytes) 
IBM_Long Floating ( 8 bytes) 
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Values for x_type and y_type Floating-Point Data Type 


CVT_CRAY_SINGLE CRAY _Floating ( 8 bytes) 


Provide a zero (0) value to the options argument to select default behavior 

or choose one or more options (status condition option, rounding options, 
"FORCE" options, CRAY and |BM options) from the tables below as the options 
argument. Specify only the options that apply to your conversion. A conflicting 
or incompatable options argument will be reported as an error (CVT_INVALID_ 
OPTION). 


Applicable 

Conversion Status Condition Option Description 

All CVT_REPORT_ALL Report all applicable status 
conditions as the default. The 
reporting of recoverable status 
conditions is disabled by default 
when this option is not used. 

Applicable 

Conversion Rounding Options Description 

All CVT_ROUND_TO_NEAREST The default rounding option for 


conversions to IEEE data types. 
This IEEE Std. 754 rounding 
mode results in the representable 
output value nearest to the 
infinitely precise result. If the two 
nearest representable values are 
equally near, the one with its least 
significant bit zero is the result. 


All CVT_BIASED_ROUNDING The default rounding option for 
conversions to non-IEEE data 
types. Performs "traditional" style 
rounding. This mode results in the 
representable output value nearest 
to the infinitely precise result. If the 
two nearest representable values are 
equally near, the result is the value 
with the largest magnitude. 


All CVT_ROUND_TO_ZERO Round the output value toward zero 
(truncate). 
All CVT_ROUND_TO_POS Round the output value toward 


positive infinity. 
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Applicable 

Conversion Rounding Options Description 

All CVT_ROUND_TO_NEG Round the output value toward 

negative infinity. 

Applicable 

Conversion "FORCE" Options Description 

All CVT_FORCE_ALL_SPECIAL_VALUES Apply all applicable "FORCE" 
options for the current conversion. 

IEEE CVT_FORCE_DENORM_TO ZERO! Force a denormalized IEEE output 
value to zero. 

IEEE CVT_FORCE_INF_TO_MAX_FLOAT! Force a positive |EEE infinite output 
value to +max_float and force a 
negative IEEE infinite output value 
to -max_float. 

IEEE or CVT_FORCE_INVALID_TO_ZERO2 Force an invalid IEEE NaN (nota 

VAX number) output value or a VAX ROP 


(reserved operand) output value to 
zero. 


1This option is valid only for conversions to |EEE output values. 
2This option is valid only for conversions to|EEE or VAX output values. 


Applicable 

Conversion Options for CRAY Format Conversion Description 

CRAY CVT_ALLOW_OVRFLW_RANGE_VALUES Allow an input/output 
exponent value > 60000 
(8). 

CRAY CVT_ALLOW_UDRFLW_RANGE_VALUES Allow an input/output 
exponent value < 20000 
(8). 

Applicable 

Conversion Option for IBM Format Conversion Description 

IBM CVT_ALLOW_UNNORMALIZED_VALUES Allow unnormalized input 


arguments. Allow an 
unnormalized output value 
for a small value which would 
normalize to zero. 


CPML-30 CPML Routines 


CPML Routines 


cvt_ftof - Convert Between Supported Floating-Point Data Types 


Returns 


The return value is a bit field containing the condition codes raised by the 
function. cvt_ftof() returns CVT_NORMAL; otherwise, it sets one or more 
of the following recoverable and unrecoverable conditions. Use the following 
condition names to determine which conditions are set: 


Condition Name 


Condition (Always reported by default) 


CVT_INVALID_INPUT_TYPE 
CVT_INVALID_OUTPUT_TYPE 
CVT_INVALID_OPTION 


Invalid input type code. 
Invalid output type code. 


Invalid option argument. 


Condition Name 


Condition (Only reported if the CVT_ 
REPORT_ALL option is selected) 


CVT_RESULT_INFINITE 
CVT_RESULT_DENORMALIZED 
CVT_RESULT_OVERFLOW_RANGE 
CVT_RESULT_UNDERFLOW_RANGE 
CVT_RESULT_UNNORMALIZED 
CVT_RESULT_INVALID 


CVT_RESULT_OVERFLOW 
CVT_RESULT_UNDERFLOW 
CVT_RESULT_INEXACT 


Conversion produced an infinite result.? 
Conversion produced a denormalized result.+ 
Conversion yielded an exponent > 60000 (8).2 
Conversion yielded an exponent < 20000 (8).2 
Conversion produced an unnormalized result.3 


Conversion result is either ROP (reserved operand), 
NaN (not a number), or closest equivalent. CRAY 
and IBM data types return 0.4 


Conversion resulted in overflow.4 
Conversion resulted in underflow.’ 


Conversion resulted in a loss of precision.* 


1For IEEE data type conversions. 
2For CRAY data type conversions. 
3For IBM data type conversions. 
4F or all data type conversions. 
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See Also 
Appendix A, Critical Floating-Point Values 


ANSI/IEEE Std 754-1985, IEEE Standard for Binary Floating-Point 
Arithmetic 


Example 


status = cvt_ftof( s&big x, CVT_BIG ENDIAN IEEE T, 
&little x, CVT_IEEE T, 0 ); 


This example converts the value pointed to by big_x, which is of type IEEE 
Big Endian T Floating, to the IEEE Little Endian T Floating data type. It 
stores the result in the location pointed to by little x. No conversion options 
are specified. 


status = cvt_ftof(&x, CVT_VAX_D, &y, CVT_IEEE T, 
(CVT_FORCE ALL SPECIAL VALUES | CVT REPORT ALL) ); 


This example converts the value pointed to by x, which is of type VAX D 
Floating, to the IEEE Little Endian T Floating data type. It stores the result 
in the location pointed to by y. Any special IEEE values that would normally 
be generated will be removed. That is, NaN and Denormalized results will be 
returned as zero and infinite results will go to + max_float. In addition, all 
recordable status conditions will be reported. 
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drem - Remainder 


Interface 
F_TYPE drem (F_TYPE x, F_TYPE y) 


Description 


drem() returns the remainder r =x-n*y, where n =rint(x/y). Additionally, if 

| n-x/y| =1/2, then n is even. The remainder is computed exactly, and | r| is 
less than or equal to | y| /2. The drem() and remainder() functions are aliases 
of each other. 


IEEE Std 754 defines drem(x,0) and drem(infinity,y) to be invalid operations 
that produce a NaN. 


Exceptions 
Exceptional Argument Routine Behavior 
y =0 Invalid argument 
x =infinity Invalid argument 
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erf - Error Functions 


Interface 
F_TYPE erf (F_TYPE x) 
F_TYPE erfc (F_TYPE x) 


Description 


erf() returns the value of the error function where erf(x) equals (2 * sqrt(pi)) 
times the area under the curve e**(-t**2) between 0 and x. 


erfc() returns (1.0-erf(x)). The erfc() function can result in an underflow as x 
gets large. 


Exceptions 


None. 
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exp - Exponential 


Interface 
F_TYPE exp (F_TYPE x) 
F_TYPE expm1 (F_TYPE x) 


Description 


exp() computes the value of the exponential function, defined as e**x, where e 
is the constant used as a base for natural logarithms. 


expm1() computes exp(x)-1 accurately, even for tiny x. 


Exceptions 
Exceptional Argument Routine Behavior 
xX >In(max_float) Overflow 
x <In(min_float) Underflow 

See Also 


Appendix A, Critical Floating-Point Values 
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fabs - Absolute Value 


Interface 
F_TYPE fabs (F_TYPE x) 


Description 


fabs( ) computes the absolute value of x. 


Exceptions 


None. 
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finite - Check for Finite Value 


Interface 
int finite (F_TYPE x) 


Description 
finite() returns the integer value 1 (true) or O (false). 
finite(x) = 1 when -infinity <x <-+Hnfinity. 


finite(x) =O when | x| =infinity or x isa NaN. 


Exceptions 


None. 
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CPML Routines 
floor - Floor 


floor - Floor 


Interface 
F_TYPE floor (F_TYPE x) 


Description 


floor() returns the largest floating-point integer value less than or equal to x. 


Exceptions 


None. 
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CPML Routines 
fmod - Modulo Remainder 


fmod - Modulo Remainder 


Interface 
F_TYPE fmod (F_TYPE x, F_TYPE y) 


Description 


fmod() computes the floating-point remainder of x modulo y. It returns the 
remainder r =x-n*y, where n =trunc(x/y). The remainder is computed exactly. 


The result has the same sign as x and a magnitude less than the magnitude of 


y. 

Exceptions 
Exceptional Argument Routine Behavior 
x =infinity Invalid argument 
y =0 Invalid argument 
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CPML Routines 
fp_class - Classifies IEEE Floating-Point Values 


fp_class - Classifies IEEE Floating-Point Values 


Interface 
int fo_class (F_TYPE x) 


Description 


These routines determine the class of IEEE floating-point values. They return 
one of the constants in the file <fp _class.h> and never cause an exception, 
even for signaling NaNs. These routines implement the recommended function 
class(x) in the appendix of the IEEE Std 754. The constants in <fp_dass.h> 
refer to the following classes of values: 


Constant Class 

FP_SNAN Signaling NaN (Not-a-Number) 
FP_QNAN Quiet NaN (Not-a-Number) 
FP_POS INF H nfinity 

FP_NEG INF -I nfinity 

FP_POS NORM Positive normalized 


FP_NEG NORM 
FP_POS DENORM 
FP_NEG DENORM 


Negative normalized 
Positive denormalized 
Negative denormalized 


FP_POS ZERO +0.0 (positive zero) 
FP_NEG ZERO -0.0 (negative zero) 
Exceptions 
None. 
See Also 


ANSI/IEEE Std 754-1985, IEEE Standard for Binary Floating-Point Arithmetic 
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CPML Routines 
frexp - Convert to Fraction and Integral Power of 2 


frexp - Convert to Fraction and Integral Power of 2 


Interface 
F_TYPE frexp (F_TYPE x, int *n) 


Description 


frexp() breaks a floating-point number into a normalized fraction and an 
integral power of 2. It stores the integer in the int object pointed to by the n 
parameter and returns the fraction part. 


Exceptions 
Exceptional Argument Routine Behavior 
| x| =infinity Invalid argument 
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CPML Routines 
hypot - Euclidean Distance 


hypot - Euclidean Distance 


Interface 
F_TYPE hypot (F_TYPE x, F_TYPE y) 


Description 


hypot() computes the length of the hypotenuse of a right triangle, where x and 
y represent the perpendicular sides of the triangle. 


hypot(x,y) is defined as the square root of (x**2 + y**2) and returns the same 
value as cabs(x,y). 


Exceptions 
Exceptional Argument Routine Behavior 
sqrt(x**2 + y**2) >max_float Overflow 

See Also 
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CPML Routines 
ilogb - Computes an Unbiased Exponent 


ilogb - Computes an Unbiased Exponent 


Interface 
int logo (F_TYPE x) 


Description 


ilogb() returns the integral part of lo (#l) as a signed integral value, for 
x iS not equal to 0, where r is the radix of the machine's floating point 


arithmetic. 
Exceptions 
Exceptional Argument Routine Behavior 
| x| =infinity INT_MAX 
x =0, NaN INT MIN 
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CPML Routines 
isnan - Check for NaN Value 


isnan - Check for NaN Value 


Interface 
intisnan (F_TYPE x) 


Description 


isnan() returns 1 (true) if x is NaN (the IEEE floating-point reserved Not-a- 
Number value) and 0 (false) otherwise. 


Exceptions 


None. 
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CPML Routines 
Idexp - Multiply by an Integral Power of 2 


Idexp - Multiply by an Integral Power of 2 


Interface 
F_TYPE ldexp (F_TYPE x, int n) 


Description 
Idexp() multiplies a floating-point number, x, by 2**n. 


Exceptions 
Exceptional Argument Routine Behavior 
| x*(2**n)| > max_float Overflow 
| x*(2**n)| <min_float Underflow 

See Also 


Appendix A, Critical Floating-Point Values 
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CPML Routines 
Igamma - Computes the Logarithm of the gamma Function 


Igamma - Computes the Logarithm of the gamma Function 


Interface 
F_TYPE lgamma (F_TYPE x) 


Description 


Igamma() returns the logarithm of the absolute value of gamma of x, or 

In(| G(x)| ), where G is the gamma function. The sign of gamma of x is 
returned in the external integer variable signgam as +1 or -1. The x parameter 
cannot be 0 or a negative integer. 


gamma() returns the natural log of the gamma function and so is functionally 
equivalent to Igamma(). Because of this, gamma() is marked TO BE 
WITHDRAWN in the X/ Open Portability Guide Revision 4 (XPG4). 


Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
x =0, -1, -2, -3, ... Invalid argument 
| x|_ >lgamma_max_float Overflow 

See Also 


Appendix A, Critical Floating-Point Values 
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CPML Routines 
log - Logarithm Functions 


log - Logarithm Functions 


Interface 
F_TYPE In (F_TYPE x) 
F_TYPE log2 (F_TYPE x) 
F_TYPE logi0 (F_TYPE x) 
F_TYPE logip (F_TYPE y) 


Description 
In() computes the natural (base e) logarithm of x. 
log2() computes the base 2 logarithm of x. 
1og10( ) computes the common (base 10) logarithm of x. 
loglp() computes In(1+y) accurately, even for tiny y. 


Exceptions 
Exceptional Argument Routine Behavior 
x <0 Invalid argument 
x =0 Overflow 
lty <0 Invalid argument 
l+ty =0 Overflow 
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CPML Routines 
logb - Radix-independent Exponent 


logb - Radix-independent Exponent 


Interface 
F_TYPE logb (F_TYPE x) 


Description 


logb() returns a signed integer converted to double-precision floating-point and 
so chosen that 1 <= | x| /2**n <2 unless x =Oor | x| =infinity. 

IEEE Std 754 defines logb(+infinity) = +infinity and logb(0) =-infinity. The 
latter is required to signal division by zero. 


Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
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CPML Routines 
modf - Return the Fractional Part and Integer Part of a Floating-Point Number 


modf - Return the Fractional Part and Integer Part of a Floating- 
Point Number 


Interface 
F_TYPE modf (F_TYPE x, F_TYPE *n) 


Description 


modf() splits a floating-point number x into a fractional part f and an integer 
part i such that | f| <1.0 and (f +i) =x. Both f and i have the same sign as 
x. modf() returns f and stores i into the location pointed to by n. 


Exceptions 


None. 
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CPML Routines 
nextafter - Next Machine Number After 


nextafter - Next Machine Number After 


Interface 
F_TYPE nextafter (F_TYPE x, F_TYPE y) 


Description 
nextafter() returns the machine-representable number next to x in the 
direction y. 
Exceptions 
Exceptional Argument Routine Behavior 
xX =max_float and y = +nfinity Overflow 
xX =-max_float and y =-infinity Overflow 
xX =min_float and y is less than or equal to 0 Underflow 
xX =-min_float and y is greater than or equal to 0 Underflow 
See Also 


ANSI/IEEE Std 754-1985, [EEE Standard for Binary Floating-Point Arithmetic 
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CPML Routines 
nint - Round to the Nearest Integer 


nint - Round to the Nearest Integer 


Interface 
F_TYPE nint (F_TYPE x) 


Description 


nint() returns the nearest integral value to x, except halfway cases are rounded 
to the integral value larger in magnitude. This function corresponds to the 
Fortran generic intrinsic function nint(). 


Exceptions 


None. 
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CPML Routines 
pow - Raise the Base to a Floating-Point Exponent 


pow - Raise the Base to a Floating-Point Exponent 


Interface 
F_TYPE pow (F_TYPE x, F_TYPE y) 


Description 


pow() raises a floating-point base x to a floating-point exponent y. The value 
of pow(x,y) is computed as e**(y In(x)) for positive x. If x is 0 or negative, see 
your language reference manual. 


Passing a NaN input value to pow() produces a NaN result for nonzero values 
of y. For pow(NaN,0), see your language reference manual. 


Exceptions 
Exceptional Argument Routine Behavior 
y In(x) >In(max_float) Overflow 
y In(x) <In(min_float) Underflow 
Fortran-Exceptional Argument Routine Behavior 
x <0 Invalid argument 
x =O and y <0 Invalid argument 
x =O and y =0 Invalid argument 
x =+infinity and y =0 Invalid argument 
xX =land | y| =infinity Invalid argument 
ANSI C-Exceptional Argument Routine Behavior 
| x| =land | y| =infinity Invalid argument 
x <O and y is not integral Invalid argument 
See Also 


Appendix A, Critical Floating-Point Values 


CPML-52 CPML Routines 


CPML Routines 
random - Random Number Generator, Uniformly Distributed 


random - Random Number Generator, Uniformly Distributed 


Interface 
F_TYPE random (int *n) 


Description 


random() is a general random number generator. The argument to the random 
function is an integer passed by reference. There are no restrictions on 

the input argument, although it should be initialized to different values on 
separate runs in order to obtain different random sequences. This function 
must be called again to obtain the next pseudorandom number. The argument 
is updated automatically. 

The result is a floating-point number that is uniformly distributed in the 
interval (0.0,1.0). 


Exceptions 
None. 
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CPML Routines 
remainder - Remainder 


remainder - Remainder 


Interface 
F_TYPE remainder (F_TYPE x, F_TYPE y) 


Description 


remainder() returns the remainder r =x-n*y, where n =rint(x/y). Additionally, 
if | n-x/y| =1/2, then n is even. Consequently, the remainder is computed 
exactly, and | r| is less than or equal to | y| /2. The drem() and remainder() 
functions are aliases of each other. 


IEEE Std 754 defines remainder(x,0) and remainder(infinity,y) to be invalid 
operations that produce a NaN. 


Exceptions 
Exceptional Argument Routine Behavior 
y =0 Invalid argument 
x =infinity Invalid argument 
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CPML Routines 
rint - Return the Nearest Integral Value 


rint - Return the Nearest Integral Value 


Interface 
F_TYPE rint (F_TYPE x) 


Description 


rint() rounds x to an integral value according to the current IEEE rounding 
direction specified by the user. 


Exceptions 


None. 


CPML Routines CPML-55 


CPML Routines 
scalb - Exponent Adjustment 


scalb - Exponent Adjustment 


Interface 
F_TYPE scalb (F_TYPE x, F_TYPE y) 


Description 
scalb() = x*(2**y) computed, for integer y. 


Exceptions 
Exceptional Argument Routine Behavior 
x*(2**y) > max_float Overflow 
x*(2**y) <min_float Underflow 

See Also 
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CPML Routines 
sin - Sine of Angle 


sin - Sine of Angle 


Interface 
F_TYPE sin (F_TYPE x) 
F_TYPE sind (F_TYPE x) 


Description 
sin() computes the sine of x, measured in radians. 


sind() computes the sine of x, measured in degrees. 


Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
(sind) | x| <(180/pi) * min_float Underflow 

See Also 


Appendix A, Critical Floating-Point Values 


CPML Routines CPML-—57 


CPML Routines 
sincos - Sine and Cosine of Angle 


sincos - Sine and Cosine of Angle 


Interface 
F_COMPLEX sincos (F_TYPE x) 
F_COMPLEX sincosd (F_TYPE x) 


Description 
sincos() computes both the sine and cosine of x, measured in radians. 
sincosd() computes both the sine and cosine of x, measured in degrees. 
sincos(x) is defined as (sin x +icos y). 


Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
(sind) | x] <(180/pi) * min_float Underflow 
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CPML Routines 
sinh - Hyperbolic Sine 


sinh - Hyperbolic Sine 


Interface 
F_TYPE sinh (F_TYPE x) 


Description 
sinh() computes the hyperbolic sine of x. 
sinh(x) is defined as (exp(x)-exp(-x))/2. 


Exceptions 
Exceptional Argument Routine Behavior 
| x] >In(2 * max_float) Overflow 

See Also 


Appendix A, Critical Floating-Point Values 


CPML Routines CPML-59 
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sinhcosh - Hyperbolic Sine and Cosine 


sinhcosh - Hyperbolic Sine and Cosine 


Interface 
F_COMPLEX sinhcosh (F_TYPE x) 


Description 
sinhcosh( ) computes both the hyperbolic sine and hyperbolic cosine of x. 


sinhcosh(x) is defined as (sinh x +icosh x). 


Exceptions 
Exceptional Argument Routine Behavior 
| x] >In(2 * max_float) Overflow 

See Also 
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CPML Routines 
sqrt - Square Root 


sqrt - Square Root 


Interface 
F_TYPE sqrt (F_TYPE x) 


Description 
sqrt() computes the rounded square root of x. 


For platforms supporting a signed zero, sqrt(-0) =0. 


Exceptions 
Exceptional Argument Routine Behavior 
x <0 Invalid argument 
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CPML Routines 
tan - Tangent of Angle 


tan - Tangent of Angle 


Interface 
F_TYPE tan (F_TYPE x) 
F_TYPE tand (F_TYPE x) 


Description 
tan() computes the tangent of x, measured in radians. 


tand() computes the tangent of x, measured in degrees. 


Exceptions 
Exceptional Argument Routine Behavior 
| x] =infinity Invalid argument 
(tand) | x] <(180/pi) * min_float Underflow 
(tand) x =(2n+1) * 90 Overflow 

See Also 
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tanh - Hyperbolic Tangent 


tanh - Hyperbolic Tangent 


Interface 
F_TYPE tanh (F_TYPE x) 


Description 
tanh() computes the hyperbolic tangent of x. 
tanh(x) is defined as (exp(x)-exp(-x))/(exp(x) + exp(-x)). 


Exceptions 


None. 
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CPML Routines 
trunc - Truncation 


trunc - Truncation 


Interface 
F_TYPE trunc (F_TYPE x) 


Description 


trunc() truncates x to an integer. 


Exceptions 


None. 
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CPML Routines 
unordered - Check for x Unordered with Respect to y 


unordered - Check for x Unordered with Respect to y 


Interface 
int unordered (F_TYPE x, F_TYPE y) 


Description 


unordered(x,y) returns the value 1 (true) if x, y, or both, area NaN and returns 
the value 0 (false) otherwise. 


Exceptions 


None. 
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A 


Critical Floating-Point Values 


Table A-1 contains the hexadecimal and decimal boundary values used in 
CPML calculations and exception checking. 


Table A-1 Hexadecimal and Decimal Boundary Values 


FFFFIFFF 
FFFFFFFFFFFF7FFF 
7F 7FF FFF 
7 EFFFFFFFFFFFFF 


Data 

Type Value for: max_float 
F Hexadecimal: 

G Hexadecimal: 

iS) Hexadecimal: 

T Hexadecimal: 

Xx Hexadecimal: 


7F FEF FFFRFRFFRFRRERRFERPPERPPFFFFE 


<x oI WOT 


Decimal: 1.701411e38 

Decimal: 8.988465674311579e307 

Decimal: 3.402823e38 

Decimal: 1.797693134862316e308 

Decimal: 1.189731495357231765085759326628007016196477e4932 


(continued on next page) 
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Critical Floating-Point Values 


Table A-1 (Cont.) Hexadecimal and Decimal Boundary Values 


00000080 
0000000000000010 
00000001 
0000000000000001 


Data 

Type Value for: min_float 
F H exadecimal: 

G Hexadecimal: 

S Hexadecimal: 

T H exadecimal: 

Xx H exadecimal: 


00000000000000000000000000000001 


Decimal: 2.9387359e-39 
Decimal: 5.562684646268003e-309 


Decimal: 4.940656458412465e-324 


F 
G 
S) Decimal: 1.4012985e-45 
T 
X 


Decimal: 6.4751751194380251109244389582276465524996e-4966 


OF 3443B0 
7B 616E 3A28B 740A6 
42B17218 
40862E 42F EFA39E F 


Data 

Type Value for: In(max_float) 
F H exadecimal: 

G H exadecimal: 

S H exadecimal: 

T H exadecimal: 

Xx H exadecimal: 


400C62E 42F E FA39E F 35793C7673007E 6 


<x AWN@A TT 


Decimal: 88.029692 

Decimal: 709.0895657128241 

Decimal: 88.7228391 

Decimal: 709.7827128933840 

Decimal: 11356.5234062941439494919310779707648912527 
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(continued on next page) 


Critical Floating-Point Values 


Table A-1 (Cont.) Hexadecimal and Decimal Boundary Values 


Data 

Type Value for: In(min_float) 

F Hexadecimal: 7218C3B1 

G Hexadecimal: 39EFFEFA2E42CO0A6 

S Hexadecimal: C2CE8ED0 

T Hexadecimal: C0874385446D71C3 

X Hexadecimal: C00C6546282207802C 89D 24D65E 96274 

F Decimal: -88.72284 

G Decimal: -709.7827128933840 

S) Decimal: -103.2789 

T Decimal: -744.4400719213813 

X Decimal: -11432.7695961557379335278266113311643138373 
Data 

Type Value for: In(2 * max_float) 

F Hexadecimal: 721843B1 

G Hexadecimal: 39E FFEFA2E 4240A6 

S Hexadecimal: 42B2D4FC 

T Hexadecimal: 408633CE 8F B9F 87E 

X Hexadecimal: 400C62E 9B B80635D81D36125B 64DA4A6 
F Decimal: 88.72284 

G Decimal: 709.7827128933840 

S) Decimal: 89.41599 

T Decimal: 710.4758600739439 

Xx Decimal: 11357.2165534747038948013483100922230678208 


(continued on next page) 
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Critical Floating-Point Values 


Table A—1 (Cont.) Hexadecimal and Decimal Boundary Values 


Data 

Type Value for: (180/pi) * min_float 

F Hexadecimal: 2EE 10365 

G Hexadecimal: C1F 81A63A5DCO006C 

S Hexadecimal: 00000039 

T Hexadecimal: 0000000000000039 

X Hexadecimal: 00000000000000000000000000000039 

F Decimal: 1.683772e-37 

G Decimal: 3.187183529933798e-307 

S) Decimal: 8.028849e-44 

T Decimal: 2.830787630910868e-322 

X Decimal: 3.71000205951917569316937757202433432154392e-4964 
Data 

Type Value for: Igamma_max_float 

F Hexadecimal: 50F 97CC6 

G Hexadecimal: F55F C5015ABD7F 67 

S Hexadecimal: 7BC650F 9 

T Hexadecimal: 7F 475ABDC501F 55F 

X Hexadecimal: 7F F171AA9917F FF BD7E A44AE 6D203DF 6 
E Decimal: 2.0594342e36 

G Decimal: 1.2812545499066958e305 

S Decimal: 2.0594342e36 

T Decimal: 1.2812545499066958e305 

X Decimal: 1.0485738685148938358098967157129705040168e4928 
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B 


CPML Entry-Point Names 


Each entry-point name in Table B-1 is unique and corresponds to data-type 
specific calculations in a CPML routine. For example, the acos function has 
five entry-point-names for the OpenVMS Alpha operating system. Because 
five floating-point data types are available, five acos routines are provided: 
math$acos_s, math$acos t, math$acos f, math$acos_g, and math$acos x. Use 
the entry-point name that corresponds to your input argument data type. 


Table B—1_ Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 
Function Data Type 
Name Required OpenVMS Alpha Tru64 UNIX Alpha 
acos S_FLOAT math$acos_s acosf 
T_FLOAT math$acos_t acos 
X_FLOAT math$acos_x 
F_FLOAT math$acos_f 
G_ FLOAT math$acos_g 
acosd S_ FLOAT math$acosd_s acosdf 
T_FLOAT math$acosd_t acosd 
X_FLOAT math$acosd_x 
F_FLOAT math$acosd_f 
G_ FLOAT math$acosd_g 


(continued on next page) 
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CPML Entry-Point Names 


Table B—1 (Cont.) Enitry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

acosh S_ FLOAT math$acosh_s acoshf 
T_FLOAT math$acosh_t acosh 
X_FLOAT math $acosh_x 
F_FLOAT math$acosh_f 
G_FLOAT math$acosh_g 

asin S_ FLOAT math$asin_s asinf 
T_FLOAT math$asin_t asin 
X_FLOAT math$asin_x 
F_FLOAT math$asin_f 
G_ FLOAT math$asin_g 

asind S_ FLOAT math$asind_s asindf 
T_FLOAT math$asind_t asind 
X_FLOAT math$asind_x 
F_FLOAT math$asind_f 
G_ FLOAT math$asind_g 

asinh S_ FLOAT math$asinh_s asi nhf 
T_FLOAT math$asinh_t asinh 
X_FLOAT math $asinh_x 
F_FLOAT math$asinh_f 
G_FLOAT math$asinh_g 
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(continued on next page) 


CPML Entry-Point Names 


Table B-1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

atan S_ FLOAT math$atan_s atanf 
T_FLOAT math$atan_t atan 
X_FLOAT math$atan_x 
F_FLOAT math$atan_f 
G_ FLOAT math$atan_g 

atan2 S_ FLOAT math$atan2_s atan2f 
T_FLOAT math$atan2_t atan2 
X_FLOAT math$atan2_x 
F_FLOAT math$atan2_f 
G_ FLOAT math$atan2_g 

atand S_ FLOAT math$atand_s atandf 
T_FLOAT math$atand_t atand 
X_FLOAT math$atand_x 
F_FLOAT math$atand_f 
G_ FLOAT math$atand_g 

atand2 S_ FLOAT math$atand2_s atand2f 
T_FLOAT math$atand2_t atand2 
X_FLOAT math$atand2_x 
F_FLOAT math$atand2_f 
G_ FLOAT math$atand2_g 


(continued on next page) 
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Table B—1 (Cont.) Enitry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

atanh S_ FLOAT math$atanh_s atanhf 
T_FLOAT math$atanh_t atanh 
X_FLOAT math$atanh_x 
F_FLOAT math$atanh_f 
G_FLOAT math$atanh_g 

cabs S_ FLOAT math$hypot_s cabsf 
T_FLOAT math$hypot_t cabs 
X_FLOAT math$hypot_x 
F_FLOAT math$hypot_f 
G_FLOAT math$hypot_g 

cort S_ FLOAT math $cbrt_s cbrtf 
T_FLOAT math$cbrt_t cort 
X_FLOAT math$cbrt_x 
F_FLOAT math$cbrt_f 
G_FLOAT math$cbrt_g 

ccos S_ FLOAT math $ccos_s ccosf 
T_FLOAT math$ccos_t ccos 
X_FLOAT math $ccos_x 
F_FLOAT math$ccos_f 
G_FLOAT math$ccos_g 
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Table B-1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

cdiv S_ FLOAT math$cdiv_s cdivf 
T_FLOAT math$cdiv_t cdiv 
X_FLOAT math$cdiv_x 
F_FLOAT math$cdiv_f 
G_ FLOAT math$cdiv_g 

ceil S_ FLOAT math$ceil_s ceilf 
T_FLOAT math$ceil_t ceil 
X_FLOAT math$ceil_ x 
F_FLOAT math$ceil_ f 
G_ FLOAT math$ceil_g 

cexp S_ FLOAT math$cexp_s cexpf 
T_FLOAT math$cexp_t cexp 
X_FLOAT math$cexp_x 
F_FLOAT math$cexp_f 
G_ FLOAT math$cexp_g 

clog S_ FLOAT math$clog_s clogf 
T_FLOAT math$clog_t clog 
X_FLOAT math$clog_x 
F_FLOAT math$clog_f 
G_ FLOAT math$clog_g 


(continued on next page) 
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Table B—1 (Cont.) Enitry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

cmul S_ FLOAT math$cmul_s cmulf 
T_FLOAT math$cmul_t cmul 
X_FLOAT math$cmul_x 
F_FLOAT math$cmul_f 
G_FLOAT math$cmul_g 

copysign S_ FLOAT math$copysign_s copysi gnf 
T_FLOAT math$copysign_t copysign 
X_FLOAT math$copysign_x 
F_FLOAT math$copysign_f 
G_FLOAT math$copysign_g 

cos S_ FLOAT math $cos_s cosf 
T_FLOAT math$cos_t cos 
X_FLOAT math$cos_x 
F_FLOAT math$cos_f 
G_FLOAT math$cos_g 

cosd S_ FLOAT math$cosd_s cosdf 
T_FLOAT math$cosd_t cosd 
X_FLOAT math$cosd_x 
F_FLOAT math$cosd_f 
G_FLOAT math$cosd_g 
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Table B-1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

cosh S_ FLOAT math$cosh_s coshf 
T_FLOAT math$cosh_t cosh 
X_FLOAT math$cosh_x 
F_FLOAT math$cosh_f 
G_ FLOAT math$cosh_g 

cot S_ FLOAT math$cot_s cotf 
T_FLOAT math$cot_t cot 
X_FLOAT math$cot_x 
F_FLOAT math$cot_f 
G_ FLOAT math$cot_g 

cotd S_ FLOAT math$cotd_s cotdf 
T_FLOAT math$cotd_t cotd 
X_FLOAT math$cotd_x 
F_FLOAT math$cotd_f 
G_ FLOAT math$cotd_g 

cpow S_ FLOAT math$cpow_s cpowf 
T_FLOAT math$cpow_t cpow 
X_FLOAT math$cpow_x 
F_FLOAT math$cpow_f 
G_ FLOAT math$cpow_g 
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Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

csin S_ FLOAT math$csin_s csinf 
T_FLOAT math$csin_t csin 
X_FLOAT math$csin_x 
F_FLOAT math $csin_f 
G_FLOAT math$csin_g 

csqrt S_ FLOAT math$csqrt_s csqrtf 
T_FLOAT math$csqrt_t csqrt 
X_FLOAT math$csqrt_x 
F_FLOAT math$csqrt_f 
G_FLOAT math$csqrt_g 

cvt_ftof All supported cvt_ftof 
types 

drem S_ FLOAT math$rem_s dremf 
T_FLOAT math$rem_t drem 
X_FLOAT math$rem_x 
F_FLOAT math$rem_f 
G_FLOAT math$rem_g 
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Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

erf S FLOAT math$erf_s erff 
T_FLOAT math$erf_t erf 
X_FLOAT math$erf_x 
F_FLOAT math$erf_f 
G_ FLOAT math$erf_g 

erfc S_ FLOAT math$erfc_s erfcf 
T_FLOAT math$erfc_t erfc 
X_FLOAT math$erfc_x 
F_FLOAT math$erfc_f 
G_ FLOAT math$erfc_g 

exp S_ FLOAT math$exp_s expf 
T_FLOAT math$exp_t exp 
X_FLOAT math$exp_x 
F_FLOAT math$exp_f 
G_ FLOAT math$exp_g 

expm1 S_ FLOAT math$expm1_s expm1f 
T_FLOAT math$expm1_t expm1 
X_FLOAT math$expm1_x 
F_FLOAT math$expm1_f 
G_ FLOAT math$expm1_g 


(continued on next page) 


CPML Entry-Point Names B-9 


CPML Entry-Point Names 


Table B—1 (Cont.) Enitry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

fabs S_ FLOAT math$fabs_s fabsf 
T_FLOAT math$fabs _t fabs 
X_FLOAT math$fabs_x 
F_FLOAT math$fabs_f 
G_ FLOAT math$fabs_g 

finite S_ FLOAT math$finite_s finitef 
T_FLOAT math$finite_t finite 
X_FLOAT math$finite_x 
F_FLOAT math$finite_f 
G_FLOAT math$finite_g 

floor S_ FLOAT math$floor_s floorf 
T_FLOAT math$floor_t floor 
X_FLOAT math$floor_x 
F_FLOAT math$floor_f 
G_FLOAT math$floor_g 

fmod S_ FLOAT math$mod_s fmodf 
T_FLOAT math$mod_t fmod 
X_FLOAT math$mod_x 
F_FLOAT math$mod_f 
G_FLOAT math$mod_g 
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Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

fp_class S_ FLOAT math$fp_class_s fp_classf 
T_FLOAT math$fp_class t fp_class 
X_FLOAT math$fp_class_x 
F_FLOAT math$fp_class_f 
G_ FLOAT math$fp_class_g 

frexp S_ FLOAT math$frexp_s frexpf 
T_FLOAT math$frexp_t frexp 
X_FLOAT math$frexp_x 
F_FLOAT math$frexp_f 
G_ FLOAT math$frexp_g 

hypot S_ FLOAT math$hypot_s hypotf 
T_FLOAT math$hypot_t hypot 
X_FLOAT math$hypot_x 
F_FLOAT math$hypot_f 
G_ FLOAT math$hypot_g 

ilogb S_ FLOAT math$ilogb_s ilogbf 
T_FLOAT math$ilogb_t ilogb 
X_FLOAT math$ilogb_x 
F_FLOAT math$ilogb_f 
G_FLOAT math$ilogb_g 
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Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

isnan S_ FLOAT math$isnan_s isnanf 
T_FLOAT math$isnan_t isnan 
X_FLOAT math$isnan_x 
F_FLOAT math$isnan_f 
G_FLOAT math$isnan_g 

jo S_FLOAT math$j0_s jOf 
T_FLOAT math$j0_t jo 
X_FLOAT math$j0_x 
F_FLOAT math$jO_f 
G_FLOAT math$jO_g 

jl S_FLOAT math$j1_s jif 
T_FLOAT math$j1_t jl 
X_FLOAT math$j1_x 
F_FLOAT math$j1_f 
G_FLOAT math$j1_g 

jn S_ FLOAT math$jn_s jnf 
T_FLOAT mathg$jn_t jn 
X_FLOAT math$jn_x 
F_FLOAT math$jn_f 
G_FLOAT math$jn_g 


B-12 CPML Entry-Point Names 


(continued on next page) 


CPML Entry-Point Names 


Table B-1 (Cont.) Entry-Point Names for CPML Platforms 


Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

Idexp S_ FLOAT math$ldexp_s I|dexpf 
T_FLOAT math$ldexp_t Idexp 
X_FLOAT math$ldexp_x 
F_FLOAT math$ldexp_f 
G_ FLOAT math$ldexp_g 

Igamma S_ FLOAT math$lgamma_s Igammaf 
T_FLOAT math$lgamma_t Igamma 
X_FLOAT math$lgamma_x 
F_FLOAT math$lgamma_f 
G_ FLOAT math$lgamma_g 

In S_ FLOAT math$In_s loof 
T_FLOAT math$In_t log 
X_FLOAT math$In_x 
F_FLOAT math$In_f 
G_ FLOAT math$In_g 

log2 S_ FLOAT math$log2_s log2f 
T_FLOAT math$log2_t log2 
X_FLOAT math$log2_x 
F_FLOAT math$log2_f 
G_ FLOAT math$log2_g 
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Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

loglO S_ FLOAT math$log10_s loglOf 
T_FLOAT math$log10_t log10 
X_FLOAT math$log10_x 
F_FLOAT math$log10_f 
G_FLOAT math$logl0_g 

loglp S_ FLOAT math$loglp_s loglpf 
T_FLOAT math$loglp_t loglp 
X_FLOAT math$loglp_x 
F_FLOAT math$loglp_f 
G_FLOAT math$loglp_g 

logb S_ FLOAT math$logb_s logbf 
T_FLOAT math$logb_t logb 
X_FLOAT math$logb_x 
F_FLOAT math$logb_f 
G_FLOAT math$logb_g 

modf S_ FLOAT math$modf_s modff 
T_FLOAT math$modf_t modf 
X_FLOAT math $modf_x 
F_FLOAT math$modf_f 
G_FLOAT math$modf_g 
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Entry-Point Names 


Generic 
Function Data Type 
Name Required OpenVMS Alpha Tru64 UNIX Alpha 
nextafter S_ FLOAT math$nextafter_s nextafterf 
T_FLOAT math$nextafter_t nextafter 
X_FLOAT math$nextafter_x 
F_FLOAT math$nextafter_f 
G_ FLOAT math$nextafter_g 
nint S_ FLOAT math$nint_s nintf 
T_FLOAT math$nint_t nint 
X_FLOAT math$nint_x 
F_FLOAT math$nint_f 
G_ FLOAT math$nint_g 
pow S_ FLOAT math$pow_ss powf 
T_FLOAT math$pow_tt pow 
X_FLOAT math$pow_xx 
F_FLOAT math$pow_ff 
G_ FLOAT math$pow_gg 
random S_ FLOAT math$random_|_s 
T_FLOAT 
X_FLOAT 
F_FLOAT math$random_|_f 
G FLOAT 
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Entry-Point Names 


Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

remainder S_ FLOAT math$rem_s remainderf 
T_FLOAT math$rem_t remainder 
X_FLOAT math$rem_x 
F_FLOAT math$rem_f 
G_FLOAT math$rem_g 

rint S_ FLOAT math$rint_s rintf 
T_FLOAT math$rint_t rint 
X_FLOAT math $rint_x 
F_FLOAT math$rint_f 
G_FLOAT math$rint_g 

scalb S_ FLOAT math$scalb_s scalbf 
T_FLOAT math$scalb_t scalb 
X_FLOAT math$scalb_x 
F_FLOAT math$scalb_f 
G_ FLOAT math$scalb_g 

sin S_ FLOAT math$sin_s sinf 
T_FLOAT math$sin_t sin 
X_FLOAT math$sin_x 
F_FLOAT math$sin_f 
G_FLOAT math$sin_g 
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Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

sincos S_ FLOAT math$sincos_s sincosf 
T_FLOAT math$sincos_t sincos 
X_FLOAT math$sincos_x 
F_FLOAT math$sincos_f 
G_ FLOAT math$sincos_g 

sincosd S_ FLOAT math$sincosd_s sincosdf 
T_FLOAT math$sincosd_t sincosd 
X_FLOAT math$sincosd_x 
F_FLOAT math$sincosd_f 
G_ FLOAT math$sincosd_g 

sind S_ FLOAT math$sind_s sindf 
T_FLOAT math$sind_t sind 
X_FLOAT math$sind_x 
F_FLOAT math$sind_f 
G_ FLOAT math$sind_g 

sinh S_ FLOAT math$sinh_s sinhf 
T_FLOAT math$sinh_t sinh 
X_FLOAT math$sinh_x 
F_FLOAT math$sinh_f 
G_ FLOAT math$sinh_g 
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Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

sinhcosh S_ FLOAT math$sinhcosh_s sinhcoshf 
T_FLOAT math$sinhcosh_t sinhcosh 
X_FLOAT math$sinhcosh_x 
F_FLOAT math $sinhcosh_f 
G_FLOAT math$sinhcosh_g 

sqrt S_ FLOAT math$sqrt_s sqrtf 
T_FLOAT math$sqrt_t sqrt 
X_FLOAT math$sqrt_x 
F_FLOAT math$saqrt_f 
G_FLOAT math$saqrt_g 

tan S_ FLOAT math$tan_s tanf 
T_FLOAT math$tan_t tan 
X_FLOAT math$tan_x 
F_FLOAT math$tan_f 
G_FLOAT math$tan_g 

tand S_ FLOAT math$tand_s tandf 
T_FLOAT math$tand_t tand 
X_FLOAT math$tand_x 
F_FLOAT math$tand_f 
G_FLOAT math$tand_g 
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Generic 

Function Data Type 

Name Required OpenVMS Alpha Tru64 UNIX Alpha 

tanh S_ FLOAT math$tanh_s tanhf 
T_FLOAT math$tanh_t tanh 
X_FLOAT math$tanh_x 
F_FLOAT math$tanh_f 
G_ FLOAT math$tanh_g 

trunc S_ FLOAT math$trunc_s truncf 
T_FLOAT math$trunc_t trunc 
X_FLOAT math$trunc_x 
F_FLOAT math$trunc_f 
G_ FLOAT math$trunc_g 

unordered S_ FLOAT math$unordered_s unorder edf 
T_FLOAT math$unordered_t unordered 
X_FLOAT math$unordered_x 
F_FLOAT math$unordered_f 
G_ FLOAT math$unordered_g 

yO S_ FLOAT math$y0_s yOf 
T_FLOAT math$y0_t yO 
X_FLOAT math$y0_x 
F_FLOAT math$y0 f 
G_ FLOAT math$y0_g 
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Entry-Point Names 


Generic 
Function Data Type 
Name Required OpenVMS Alpha Tru64 UNIX Alpha 
yl S FLOAT math$y1_s ylf 
T_FLOAT math$y1 t yl 
X_FLOAT math$y1_x 
F_FLOAT math$y1_f 
G_FLOAT math$yl1_g 
yn S_ FLOAT math$yn_s ynf 
T_FLOAT math$yn_t yn 
X_FLOAT math$yn_x 
F_FLOAT math$yn_f 
G_FLOAT math$yn_g 
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Glossary 


This glossary defines mathematical terms and symbolic names used in this 
manual. 


complex number 
See F_ COMPLEX. 


denormalized number 
A floating-point number with a value very close to zero. 


domain error 

An exception condition resulting from passing an argument whose value is 
outside the range of permissible values. 

exceptional argument 

Any argument value passed to a routine that does not return a meaningful 
result, or an argument defined differently for different platforms. 
F_COMPLEX 

A complex number identifier. F_COMPLEX indicates that a given routine 
returns two different values of the same floating-point data type. See 
Table 1-2 for more information. 

F_TYPE 

A floating-point number identifier. F_TYPE is used when it is not necessary to 
distinguish between the floating types. See Table 1-1 for more information. 
floating-point number 

See F_TYPE. 


Glossary-1 


HUGE_RESULT 
For VAX data types, HUGE_RESULT =max_float. 


For IEEE data types, HUGE RESULT =infinity. 


INV_RESULT 
For VAX data types, INV_RESULT =0. 
For IEEE data types, INV RESULT =NaN. 


invalid argument 
See domain error. 


max_float 


The largest finite number representable in the floating-point data types. See 
Appendix A for more information on max_float values. 


min_float 


The smallest positive normalized nonzero number representable in the 
floating-point data types. See Appendix A for more information on min_float 
values. 


NaN 


A floating-point value that is said to be “not a number” and contains an 
indeterminate quantity. 


overflow 


An exception condition caused by passing a floating-point value that is larger 
than the highest valid floating-point value. See max_float for additional 
information. 


range error 


An exception condition that occurs when a mathematically valid argument 
results in a function value that exceeds the range of representable values for 
floating-point data types. 


underflow 


An exception condition caused by passing a floating-point value that is 
lower than the lowest valid floating-point value. See min_float for additional 
information. 


Glossary—2 


A 


32-bit IEEE single-precision complex 
number, 1-4 
32-bit IEEE single precision number, 1-3 
Absolute value, CPML-36 
32-bit VAX single-precision complex number, 
1-4 
acos routine, CPML-5, B-1 
acosd routine, CPML-5, B-1 
acosdf routine, B-1 
acosf routine, B-1 
acosh routine, CPML-6, B-2 
acoshf routine, B-2 
32-bit VAX single-precision number, 1-3 
64-bit IEEE double-precision complex 
number, 1-4 
64-bit IEEE double-precision number, 1-3 
64-bit VAX double-precision complex number, 
1-4 
64-bit VAX double-precision number, 1-3 
128-bit IEEE extended-precision complex 
number, 1-4 
128-bit IEEE extended-precision number, 
1-3 
Arc cosine of angle, CPML-5 
hyperbolic, CPML-6 
Arc sine of angle, CPML-7 
Arc tangent of angle 
hyperbolic, CPML-11 
with one argument, CPML-9 
with two arguments, CPML-10 


Index 


asin routine, CPML-7, B-2 
asind routine, CPML-7, B-2 
asindf routine, B-2 

asinf routine, B-2 

asinh routine, CPML-8, B-2 
asinhf routine, B-2 

atan routine, CPML-9, B-3 
atan2 routine, 2-3, CPML-10, B-3 
atan2f routine, B-3 

atand routine, CPML-9, B-3 
atand2 routine, CPML-10, B-3 
atand2f routine, B-3 

atandf routine, B-3 

atanf routine, B-3 

atanh routine, CPML-11, B-4 
atanhf routine, B-4 


Base 10 logarithm, CPML-47 
Base 2 logarithm, CPML-47 
Bessel functions, CPML-12 
bessel routine, CPML-12 


Cc 


cabs routine, CPML-13, B-4 

cabsf routine, B-4 

cbrt routine, CPML-14, B-4 

cbortf routine, B-4 

ccos routine, CPML-15, B-4 

ccosf routine, B-4 

cdiv routine, CPML-4, CPML-16, B-5 
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cdivf routine, B-5 

ceil routine, CPML-17, B-5 

ceilf routine, B-5 

Ceiling, CPML-17 

cexp routine, CPML-18, B-5 

cexpf routine, B-5 

clog routine, CPML-19, B-5 

clogf routine, B-5 

cmul routine, CPML-20, B-6 

cmulf routine, B-6 

Common logarithm, CPML-47 

Complex absolute value, CPML-13 

Complex data types, 1-3 

Complex division, CPML-16 

Complex exponential, CPML-18 

Complex floating-point data types, 1-4 

Complex functions, 1-4 

Complex multiplication, CPML-20 

Complex natural logarithm, CPML-19 

Complex numbers, 1-3 
absolute value, CPML-13 
cosine of angle, CPML-15 
division, CPML-16 
exponential of, CPML-18 
exponentiation of, CPML-25 
multiplication, CPML-20 
natural logarithm of, CPML-19 
sine of, CPML-26 
square root of, CPML-27 

Complex power, CPML-25 

Complex square root, CPML-27 

Converting floating-point data types, 

CPML-28 

copysign routine, CPML-21, B-6 

copysignf routine, B-6 

cos routine, CPML-22, B-6 

cosd routine, CPML-22, B-6 

cosdf routine, B-6 

cosf routine, B-6 

cosh routine, CPML-23, B-7 

coshf routine, B-7 

Cosine and sine of angle, CPML-58 

Cosine of angle, CPML-22 
hyperbolic, CPML-23, CPML-60 
of a complex number, CPML-15 
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cot routine, CPML-24, B-7 
Cotangent of angle, CPML-24 
cotd routine, CPML-24, B-7 
cotdf routine, B-7 

cotf routine, B-7 

cpow routine, CPML-25, B-7 
cpowf routine, B-7 

csin routine, CPML-26, B-8 
csinf routine, B-8 

csqrt routine, CPML-27, B-8 
csqrtf routine, B-8 

Cube root, CPML-14 

cvt_ftof routine, CPML-28, B-8 


D 


Data types 
conversion of floating-point, CPML-28 
INPUT_ARG TYPE, 2-2 
RETURN_TYPE, 2-1 
types of, 1-2 

Decimal boundary values, A-1 

Denormalized numbers, 1-6 

Domain errors, 1-5 

drem routine, CPML-33, B-8 

dremf routine, B-8 

D_FLOAT data type, conversion to and from, 

CPML-28 


E 


Entry points, 2-2 
Entry-point names, B-1 

erf routine, CPML-34, B-9 
erfc routine, CPML-34, B-9 
erfcf routine, B-9 

erff routine, B-9 

Error functions, CPML-34 
Error handling, 2-2 


Errors 
domain, 1-5 
range, 1-5 


Euclidean distance, CPML-42 


Exception conditions, 1-5, 2-2 

Exception handler, 1-2 

Exceptional arguments, 1-4, 1-5, 2-3 

exp routine, CPML-35, B-9 

expf routine, B-9 

expm1 routine, CPML-35, B-9 

expm1f routine, B-9 

Exponent adjustment, CPML-56 

Exponential, CPML-35 

Exponential of a complex number, 
CPML-18 


F 


fabs routine, CPML-36, B-10 

fabsf routine, B-10 

finite routine, CPML-37, B-10 

Finite value, checks for, CPML-37 

finitef routine, B-10 

Floating-point complex data types 
for OpenVMS Alpha systems, 1-4 
for Tru64 UNIX Alpha systems, 1-4 
IEEE, 1-4 
types of, 1-4 
VAX, 1-4 

Floating-point data types 
complex, 1-3, 1-4 
conversion to other types, CPML-28 
for OpenVMS Alpha systems, 1-3 
for Tru64 UNIX Alpha systems, 1-3 
IEEE, 1-3 
types of, 1-3 
VAX, 1-3 

Floating-point number conversion 
Big Endian_IEEE_S Floating, CPML-28 
Big Endian_|EEE_T Floating, CPML-28 
Big Endian_|IEEE_X_Floating, CPML-28 
CRAY_Floating, CPML-28 
D_Floating, CPML-28 
F_Floating, CPML-28 
G Floating, CPML-28 
H_ Floating, CPML-28 
IBM_Long Floating, CPML-28 
IBM_Short_Floating, CPML-28 
IEEE_S Floating, CPML-28 


Floating-point number conversion (cont’d) 
|EEE_T Floating, CPML-28 
IEEE_X Floating, CPML-28 

floor routine, CPML-38, B-10 

floorf routine, B-10 

fmod routine, CPML-39, B-10 

fmodf routine, B-10 

fp_class routine, CPML-40, B-11 

fp_classf routine, B-11 

frexp routine, CPML-41, B-11 

frexpf routine, B-11 

F_ COMPLEX datatype, 1-3 

F_FLOAT datatype, 1-3 
conversion to and from, CPML-28 

F_FLOAT_COMPLEX data type, 1-4 


G 


gamma routine, CPML-46 
Generic interface names, 2-1 
G_FLOAT datatype, 1-3 

conversion to and from, CPML-28 
G_FLOAT_COMPLEX data type, 1-4 


H 


Hexadecimal boundary values, A-1 

Hyperbolic arc cosine of angle, CPML-6 

Hyperbolic arc sine of angle, CPML-8 

Hyperbolic arc tangent of angle, CPML-11 

Hyperbolic cosine of angle, CPML-23 

Hyperbolic sine, CPML-59 

Hyperbolic sine and cosine, CPML-60 

Hyperbolic tangent, CPML-63 

hypot routine, CPML-42, B-11 

hypotf routine, B-11 

H_FLOAT data type, conversion to and from, 
CPML-28 


IEEE Big Endian, CPML-28 
|EEE considerations, 1-6 


Index—3 


IEEE double-precision complex data type, 
1-4 

IEEE double-precision data type, 1-3 

IEEE extended-precision complex data type, 
1-4 

IEEE extended-precision data type, 1-3 

IEEE floating-point data types, conversion to 
other data types, CPML-28 

IEEE floating-point values, identifying the 
dass of, CPML-40 

IEEE Little Endian, CPML-28 

IEEE single-precision complex data type, 
1-4 

IEEE single-precision data type, 1-3 

ilogb routine, CPML-43, B-11 

ilogbf routine, B-11 

Input arguments, 2-2 

Integer data type, 1-2 

Integers, rounding, CPML-51, CPML-55 

Integral power of 2, converting to, CPML-41 

Invalid argument exception condition, 1-6 

Invalid arguments, 1-5 

isnan routine, CPML-44, B-12 

isnanf routine, B-12 


J 


jO routine, CPML-12, B-12 
jOf routine, B-12 
jl routine, CPML-12, B-12 
jlf routine, B-12 
jn routine, CPML-12, B-12 
jnf routine, B-12 


L 


Language-specific routine behavior, 1-2 
Idexp routine, CPML-45, B-13 

Idexpf routine, B-13 

Igamma routine, CPML-46, B-13 
Igammaf routine, B-13 
Igamma_max_float boundary value, A-4 
In routine, CPML-47, B-13 


Index—4 


log routine, CPML-47, B-13 
logl0 routine, CPML-47, B-14 
logl0f routine, B-14 
loglp routine, CPML-47, B-14 
loglpf routine, B-14 
log2 routine, CPML-47, B-13 
log2f routine, B-13 
Logarithm 
base 10, CPML-47 
base 2, CPML-47 
common, CPML-47 
complex number of a, CPML-19 
gamma function, CPML-46 
ilogb, compute unbiased exponent, 
CPML-43 
Igamma function, CPML-46 
logb, convert to double-precision 
floating-point, CPML-48 
natural, CPML-47 
logb routine, CPML-48, B-14 
logbf routine, B-14 
logf routine, B-13 


Machine numbers, CPML-50 
math$acosd_f, 
math$acosd_g, 
math$acosd_s, 
math$acosd t, 
math$acosd_x, 
math$acosh_f, 
math$acosh_g, 
math$acosh_s, 
math$acosh_t, 
math$acosh_x, 
math$acos f, B-1 
math$acos_g, B-1 
math$acos_s, B-1 

B-1 

B-1 
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B- 
B- 
B- 
B- 
B- 
B- 
B- 
B- 
B- 
B- 
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math$acos t, 
math$acos_x, 
math$asind_f, 


B-2 
math$asind_g, B-2 


math$asind_s, B-2 
math$asind_t, B-2 
math$asind_x, B-2 
math$asinh_f, B-2 
math$asinh_g, B-2 
math$asinh_s, B-2 
math$asinh_t, B-2 
math$asinh_x, B-2 


math$asin_f, B-2 
math$asin_g, B- 
math$asin_s, 
math$asin_t, 
math$asin_x, B- 
math$atan2_f, B-3 
math$atan2_g, B 
math$gatan2_s, B- 
math$atan2 t, B 
math$gatan2_x, B- 
math$atand2 f, B- 
math$atand2_g, 

math$atand2_s, 
math$atand2 t, 
math$atand2_x, 
math$atand_f, B- 
math$atand_g, B 
math$atand_s, B- 
mathgatand_t, B 
math$atand_x, 
math$atanh_f, 
mathgatanh_g, B 
mathgatanh_s, B- 
math$atanh_t, B 
math$atanh_x, 
mathgatan_f, 
math$atan_g, B 

math$atan_s, B-3 
mathgatan_t, B 

math$atan_x, 
math$cbrt_f, B-4 
math$cbrt_g, B 
math$cbrt_s, B 
math$cbrt_t, B- 
math$cbrt_x, B 
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math$ccos_f, 
math$ccos_g, 
math$ccos s, 
math$ccos t, 
math$ccos_x, 
math$cdiv_f, 
math$cdiv_g, 
math$cdiv_s, 
math$cdiv_t, 
math$cdiv_x, 
math$ceil_f, 
math$ceil_g, 
math$ceil_s, 
mathé$ceil_t, 
math$ceil_ x, 
math$cexp f, 
math$cexp _g, 
math$cexp s, 
math$cexp t, 
math$cexp_x, 
math$clog_f, 
math$clog_g, 
math$clog_s, 
math$clog_t, 
math$clog_x, 
math$cmul_f, 
math$cmul_g, 
math$cmul_s, 
math$cmul_t, 
math$cmul_x, 
math$copysign_ 
math$copysign_ 
math$copysign_ 
math$copysign_ 
math$copysign__ 
math$cosd_f, 
math$cosd_g, 
math$cosd_s, 
math$cosd_t, 
math$cosd_x, 
math$cosh_f, 
math$cosh_g, 
math$cosh_s, 
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math$cosh_t, B-7 
math$cosh_x, B-7 
math$cos f, B-6 

6 
math$cos_s, B-6 
math$cos t, B-6 

6 
math$cotd_f, B-7 
math$cotd_g, 
math$cotd_s, 
math$cotd t, 


math$cotd_x, B- 
math$cot_f, B-7 


math$cot_g, B-7 
math$cot_s, B-7 
math$cot_t, B-7 
math$cot_x, B-7 


math$cpow_f, B-7 
math$cpow_g, 
math$cpow_s, 
math$cpow t, 
math$cpow_x, 
math$csin_f, 
math$csin_g, 
math$csin_s, 
math$csin_t, 
math$csin_x, 
math$csaqrt_f, 
math$csqrt_g, 
math$csqrt_s, 
math$csqrt_t, 
math$csqrt_x, 
math$erfc_f, 
math$erfc_g, 
math$erfc_s, 
math$erfc t, 
math$erfc_x, 
math$erf_f, 
math$erf_g, 
math$erf_s, 
math$erf_t, 
math$erf_x, 
math$expm1_ 


18] 
1 ! 


w 
wPPOw 
wm mmm © 


mnUwom 
co 0 0 0 


ri 
o 


DUaw 
oVo0~”o 


i 
o 


DVaw 
i} I 
o 
pmoVrVvno 


hh 
1 
ie} 


’ 


Index—6 


math$expml1_g, B-9 
math$expm1_s, B-9 
math$expm1_t, B-9 
math$expm1_x, B-9 
math$exp f, B-9 


math$exp_g, B-9 
math$exp_s, B-9 
math$exp_t, B-9 


math$exp_x, B-9 
math$fabs f, B-10 
math$fabs_g, B-10 
math$fabs_s, B-10 
math$fabs_t, B-10 
math$fabs_x, B-10 
math$finite_f, B-10 
math$finite_g, B-10 
math$finite_s, B-10 
math$finite_t, B-10 
math$finite_x, B-10 
math$floor_f, B-10 
math$floor_g, B-10 
math$floor_s, B-10 
math$floor_t, B-10 
math$floor_x, B-10 
math$fp_class f, B-11 
math$fp_class_g, B-11 
math$fp_class_s, B-11 
math$fp_class t, B-11 
math$fp_class x, B-11 
math$frexp_f, B-11 
math$frexp_g, B-11 
math$frexp_s, B-11 
math$frexp_t, B-11 
math$frexp_x, B-11 
math$hypot_f, B-4, B-11 
math$hypot_g, B-4, B-11 
math$hypot_s, 4, B-11 
math$hypot_t, 4, B-11 
math$hypot_x, B-4, B-11 
math$ilogb f, B-11 
math$ilogb_g, B-11 
math$ilogb_s, B-11 
math$ilogb t, B-11 


B- 
B- 


math$ilogb_x, B-11 
math$jO_f, B-12 
math$jO_g, B-12 
math$j0_s, B-12 
mathg$jO t, B-12 
math$jO_x, B-12 
math$j1_f, B-12 
math$j1_g, B-12 
math$j1_s, B-12 
math$j1_t, B-12 
math$jn_f, B-12 
math$jn_g, B-12 
math$jn_s, B-12 
math$jn_t, B-12 
math$jn_x, B-12 
math$ldexp_f, B-13 
math$ldexp_g, B-13 
math$ldexp_s, B-13 
math$ldexp_t, B-13 
math$ldexp_x, B-13 
math$lgamma_f, B-13 
math$lgamma_g, B-13 
math$lgamma_s, B-13 
math$lgamma_t, B-13 
math$lgamma_x, B-13 
math$In_f, B-13 
math$In_g, B-13 
math$In_s, B-13 
math$In_t, B-13 
math$In_x, B-13 
math$logl0_f, B-14 
math$logl0_g, B-14 
math$logl0_s, B-14 
math$log10_t, B-14 
math$logl0_x, B-14 
math$loglp f, B-14 
math$loglp_g, B-14 
math$loglp_s, B-14 
math$loglp t, B-14 
math$loglp_x, B-14 
math$log2_f, B-13 
math$log2_g, B-13 
math$log2_s, B-13 


math$log2_t, B-13 
math$log2_x, B-13 
math$logb_ f, B-14 
math$logb_g, B-14 
math$logb_s, B-14 
math$logb_t, B-14 
math$logb_x, B-14 
math$modf_f, B-14 
math$modf_g, B-14 
math$modf_s, B-14 
math$modf_t, B-14 
math$modf_x, B-14 
math$mod_f, B-10 
math$mod_g, B-10 
math$mod_s, B-10 
math$mod_t, B-10 
math$mod_x, B-10 
math$nextafter_f, B-15 
math$nextafter_g, B-15 
math$nextafter_s, B-15 
math$nextafter_t, B-15 
math$nextafter_x, B-15 
math$nint_f, B-15 
math$nint_g, B-15 
math$nint_s, B-15 
math$nint_t, B-15 
math$nint_x, B-15 
math$pow_ff, B-15 
math$pow_gg, B-15 
math$pow_ss, B-15 
math$pow_tt, B-15 
math$pow_xx, B-15 
math$random_| f, B-15 
math$random_|_s, B-15 


math$rem_f, B-8, B-16 
math$rem_g, B-8, B-16 
math$rem_s, B-8, B-16 
math$rem_t, B-8, B-16 
math$rem_x, B-8, B-16 
math$rint_f, B-16 
math$rint_g, B-16 
math$rint_s, B-16 
math$rint_t, B-16 
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math$rint_x, B-16 
math$scalb f, B-16 
math$scalb_g, B-16 
math$scalb_s, B-16 
math$scalb_t, B-16 
math$scalb_x, B-16 
math$sincosd_f, B-17 
math$sincosd_g, B-17 
math$sincosd_s, B-17 
math$sincosd_t, B-17 
math$sincosd_x, B-17 
math$sincos f, B-17 
math$sincos_g, B-17 
math$sincos_s, B-17 
math$sincos t, B-17 
math$sincos_x, B-17 
math$sind_f, B-17 
math$sind_g, B-17 
math$sind_s, B-17 


math$sqrt_s, 

math$sqrt_t, 

math$sqrt_x, 

math$tand_f, 
math$tand_g, 
math$tand_s, 
math$tand_t, 
math$tandc_x, 
math$tanh_f, 
math$tanh_g, 
math$tanh_s, 
math$tanh_t, 
math$tanh_x, 


B-18 
B-18 
B-18 
B-18 
B-18 
B-18 
B-18 
B-18 
B-19 
B-19 
B-19 
B-19 
B-19 


math$tan_f, B-18 
math$tan_g, B-18 
math$tan_s, B-18 
math$tan_t, B-18 
math$tan_x, B-18 


math$sind_t, B-17 


math$trunc_f, B-19 
math$trunc_g, B-19 


math$sind_x, B-17 
math$sinhcosh_f, B-18 
math$sinhcosh_g, B-18 
math$sinhcosh_s, B-18 
math$sinhcosh_t, B-18 
math$sinhcosh_x, B-18 
math$sinh_f, B-17 
math$sinh_g, B-17 


math$trunc_s, B-19 
math$trunc_t, B-19 
math$trunc_x, B-19 
math$unordered_f, B-19 
math$unordered_g, B-19 
math$unordered_s, B-19 
math$unordered_t, B-19 
math$unordered_x, B-19 


math$sinh_s, B-17 
math$sinh_t, B-17 
math$sinh_x, B-17 
math$sin_f, B-16 
math$sin_g, B-16 
math$sin_s, B-16 
math$sin_t, B-16 
math$sin_x, B-16 
math$snan_f, B-12 
math$snan_g, B-12 
math$snan_s, B-12 
math$snan_t, B-12 
math$snan_x, B-12 
math$sqrt_f, B-18 
math$sqrt_g, B-18 
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math$y0 f, B-19 
math$y0_g, B-19 
math$y0_s, B-19 
math$y0 t, B-19 
math$y0_x, B-19 
math$y1_f, B-20 
math$y1_g, B-20 
math$y1_s, B-20 
math$y1_t, B-20 
math$y1_x, B-20 
math$yn_f, B-20 
math$yn_g, B-20 
math$yn_s, B-20 
math$yn_t, B-20 
math$yn_x, B-20 


max_float boundary value, A-1 
min_float boundary value, A-2 
modf routine, CPML-49, B-14 
modff routine, B-14 

Modulo remainder, CPML-39 


N 


NaN value, checking for, CPML-44 

Natural logarithm, CPML-47 
complex number of a, CPML-19 

nextafter routine, CPML-50, B-15 

nextafterf routine, B-15 

nint routine, CPML-51, B-15 

nintf routine, B-15 

Normalized fractions, converting to, 

CPML-41 


O 


Overflow exception condition, 1-6 


P 


pow routine, CPML-52, B-15 
powf routine, B-15 


R 


random routine, CPML-53, B-15 
Range errors, 1-5 
Remainder 

drem function, CPML-33 

modulo, CPML-39 

remainder function, CPML-54 
remainder routine, CPML-54, B-16 
remainderf routine, B-16 
return type, 2-1 
Right triangle, hypotenuse of a, CPML-42 
rint routine, CPML-55, B-16 
rintf routine, B-16 
Rounding to the nearest integer, CPML-51 
Routine interface, 2-1 

examples, 2-3 


S 


scalb routine, CPML-56, B-16 
scalbf routine, B-16 
sin routine, CPML-57, B-16 
sincos routine, CPML-58, B-17 
sincosd routine, CPML-58, B-17 
sincosdf routine, B-17 
sincosf routine, B-17 
sind routine, CPML-57, B-17 
sindf routine, B-17 
Sine and cosine of angle, CPML-58 
Sine of angle, CPML-57 
hyperbolic, CPML-60 
Sine of angle of a complex number, 
CPML-26 
Sine, hyperbolic, CPML-59 
sinf routine, B-16 
sinh routine, CPML-59, B-17 
sinhcosh routine, CPML-60, B-18 
sinhcoshf routine, B-18 
sinhf routine, B-17 
sqrt routine, CPML-61, B-18 
sqrtf routine, B-18 
Square root, CPML-61 
of complex numbers, CPML-27 
Symbolic constants, 1-5 
S_FLOAT data type, 1-3 
S_FLOAT data type, conversion to and from, 
CPML-28 
S_FLOAT data type, IEEE Big Endian, 
CPML-28 
S_FLOAT data type, IEEE Little Endian, 
CPML-28 
S_FLOAT_COMPLEX data type, 1-4 


T 


tan routine, CPML-62, B-18 
tand routine, CPML-62, B-18 
tandf routine, B-18 
tanf routine, B-18 
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Tangent of angle, CPML-62 
hyperbolic, CPML-63 

tanh routine, CPML-63, B-19 

tanhf routine, B-19 

trunc routine, CPML-64, B-19 

Truncation, CPML-64 

truncf routine, B-19 

T_FLOAT data type, 1-3 

T_FLOAT data type, conversion to and from, 
CPML-28 

T_FLOAT data type, [EEE Big Endian, 
CPML-28 

T_FLOAT data type, IEEE Little Endian, 
CPML-28 

T_FLOAT_COMPLEX data type, 1-4 


U 


Underflow exception condition, 1-6 
unordered routine, CPML-65, B-19 
unorderedf routine, B-19 
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V 


VAX double-precision complex data type, 
1-4 
VAX double-precision data type, 1-3 
VAX single-precision complex data type, 1-4 
VAX single-precision data type, 1-3 


X 


X/Open Portability Guide, Version 4 
conformance to, 1-7 

X_FLOAT data type, 1-3 

X_FLOAT data type, conversion to and from, 
CPML-28 

X_FLOAT data type, IEEE Big Endian, 
CPML-28 

X_FLOAT data type, IEEE Little Endian, 
CPML-28 

X_FLOAT_COMPLExX data type, 1-4 


Y 


yO routine, CPML-12, B-19 
yOf routine, B-19 
yl routine, CPML-12, B-20 
ylf routine, B-20 
yn routine, CPML-12, B-20 
ynf routine, B-20 


